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Form C.~ Computing form for Weissinger's method. 

(Underscored numbers are sample calculations.) 


O" 

--1 


RACA TN No. 14 t6 


ON 

00 


r 

2 

3 

4 

5 

6 

7 

8 

9 

10 

II 

12 

13 

14 

15 

1/ 

n 

H- 

fn^ 

Table 

CH 

(20>c, 

or 

('O)c. 


V 

n 

-I6gi/,n 

-1 ^ 
16 fiy 
-1 

l6(4)c,orc2 

{9)x(IO) 

2Bi/n 

Table 

cm 

(12) -(II) 

2byy 

Table 

CIS 

Coef. of 
’ 

(14) -(13) 

1 

1 

0 

2.613 

-2.335 


1 

1 

-1.4192 

-.3036 

.4309 

0 

-.4309 

10.4524 

10.8833 



1 

-1.414 




2 




3.8284 


0 




2 

-1.531 




3 




0 


0 




3 

.414 




4 




.2928 


0 


2 

2 

0 

-1.414 



2 

1 




2.0720 


0 




L_ 

3.696 




2 




0 


5.6568 


P.J 







3 

L 4 







Equation 

1 

2 

3 

4 

Constant 

1 

1 

1 

_L 

6, /a 

10.8833 




Gt/o 

-2.8680 




Gb / a 

.2685 




64/a 

-.1332 





(4.) 


n 

or 

1 / 

Gn Gy 

-a 

Span station 
cos-^ 

1 

.180 

.9239 

2 

.354 

.7071 

3 

.534 

.3827 

4 

.654 

0 



Form C- - Concluded. 


NACA TN No. ik-jS 
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TECHFIC/i NOTE NO. lk'j6 


A C0I4PAEI30N OF THSEE THEORETICAL METKOIB OF CALCULATING 
SPAN LOAD DISTRIBUTION ON SWEPT WINGS 
By Nicholas H. Van Dcm and John DeYoung 


SUMMARY 

Three methods for calculating span load distribution, those 
developed by V.M. Fallcner, Mutterperl, and J. Weissinger, have 
been applied to five swept wings. The angles of sweep ranged from 
—45° to +45°. These methods were examined to establish their 
relative a.cciuracy and ease of application. Experimentally determined 
loadings were used as a basis for judgin,3 accuracy. For the 
convenience of the i-eaders the computing forms and all information 
requisite to thoir application are included in appendixes. 

From the analysis it was found that the Weissinger method 
would be best suited to an over— all studjr of the effects of plan 
form on the span loading and associated characteristics of wings. 

The method gave good, but not best, accirracy and involved by far 
the least computing effort. The Falkner method gave the best 
accuracy out at a considerable expense in computing effort and 
hence appeared to be most useful for a detailed study of a specific 
wing. Tne Muttexperl method offered no advantages in accuracy or 
facility over either of the other methods. 


INTRODUCTION 

In an effort to reach higher fli^t speeds, designers are 
tixming to widely diversified types of plan fomis the aerodynamic 
characteristics of which are as yet unlmoTO. Since the multiplicity 
of such designs precludes an experimental investigation of each, 
considerable attention has been directed toward means of obtaining 
these characteristics theoretically. Usually the basis for such theo- 
retical investigations is span loading. While the precise 'span load- 
ing itself may not be considered of major importance, it is believed 
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that anj^ method giving roasonahly accm'’ate predictions of span 
loading wotxld be amenable to simple extensions which would give 
reasonably accurate values of such characteristics as lift— curve 
slope, spanwise center of pressure position, downwash at arbitrary 
locations, and rolling moments due to sideslip or rolling. 

A number of methods have been developed for predicting the 
span loading of swept wings of arbitrary taper and aspect ratio, 
but very few attempts have been ma.de to eomparo, for several methods, 
predicted and experimentally measured loadings on identical wings. 

The investiga.tlon reported herein, was undortakf;n to provide such a 
comparison of predicted a.nd measured span loadings. The theoretical 
methods have been evaluated in terras of rolacive acc\iracy, manner 
and consistency cf error, a.nd tediousneso of application. 

The methods developed by V.M. Falknor (reference l), vftn. Mutterperl 
(reference 2), and J. Weissingsr (reference 5) have been applied to 
five wings produced by sweeping the vring panels of a,n aiiplane throu^ 
a range of to +45'^. The spa.n load distributions so calculated 

have been comparcci with those obta.ined experimentally . In addition, 
the lift— curve slopes and spanwise center— of— pressure position 
predicted for each wing by the scvera.l methods have been compared 
with those vai.ues obtained experimentally. 

Throughout the calculations a chock was made of the tirnie 
required for each method and for the various parts of each method. 

From these observations a. compax'ison was made of the relative 
tediousness of oa.ch method, and Indications were obtained as to 
vfbich parts might- be rendered less diffici;.lt and time consuming. 

Finally, in order to ena.ble immediate application of the 
methods all necess.«j.’y tables, computation forms, and step-by— step 
computation instiuictions for each are included in the appendixes. 

It is believed that with these aids a computing staff could under- 
take the compxtto.tion of swept— wing chaa'-acteristics with little- 
additional supox'vision. In addition, for the convenience of the 
reader, there are included in the appendixes any ma.thematlcal 
derivations or developments not immediately obtoinable from the 
references. 
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SYMBOLS • 

General 

wing area, square- feet 

effective^ wing span, feet 

effective aspect ratio (b^/s) 

semispan (b/2), fe§t 

wing chord, feet 

root chord, feet 

average chord (S/b), feet 

taper ra.tio, tip chord divided b;' root chord (c+/co) 

sweep angle of quarter— chord line positive for sweepback, 
degrees 

geometri.c angle of attack of wing measured from an.gle for 
zero liftj degrees 

geometric angle of attack of wing root section, degrees 

local geometric angle of attack, degrees 

longitudinal coordinate of dorawash point positive 
forward, feet 

lateral coordinate of downwash point positive to 

lu.glit, feet ‘ 

dimensionless lateral cooi-dinata of doTOwash 
point (y/s) 


^In all instances except tlie unswept wing, the actual tip chord va,s 
not parallel to the wind stream. An effective tip chord that 
was parallel was therefore s.ss'umed swell that the ifing area 
remained constant. The effective span is the span to this 
eflectlve tip. 
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longitudinal coordinate of vortex olement positive 
forwai’d, feot 

lateral coordinate of vertex element positive to 
right, feet 

dimensionless lateral coordinate of vortex element (y/s) 
chord a.t spanwise station q, feet 
density of air, slugs per cuhic foot 
air— stream velocity, feet per second 

air— stream dynamic pressure (|pV 2 ), pounds per square foot 

lift, pounds 

lift coefficient (h/qS' 

section lift coefficient 

vorticity, feet per second 

circulation, feet squared per second 

spanvisG center of pressure position 

differential pressure between upper and lower surfaces 
of wing, pounds per square foot 

static pressure, poimds per square foot 

froG— stx'eara static pressm-e, pounds per sqxiaro foot 

pressure coeff icient [(p— P q) / q ] 

induced vortical velocity or downwash, feot .per second 

downwash angle, the ratio of dovmwash to free- stream 
velocity 

spanwise position in circular coordinates (cos 
or cos~^) 
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Syabols Pertaining to the Falkncr Method 

semispan of horsoshoo vortex (s/2C), foot 

dimonG ionless longitudinal coordf.nato (x/jv) 

diacnsionlcss longitudi.nal coordinate of control point 
relative to vortex 

longitudinal coordinate referred to 0.5c lino, foot 
circular longitudinal coordinate (cos~“ 2x' ,'c) 
dimoiosionloss lateral coordinato (yZ/v) 

dimons ionic S 3 lateral coordinate of control poinu rolaui ve- 
to vortex 

vinknoi-nes in distributi on series 

number of vortices in cherdvise direction 

desirgnates vhich of M vortices in chord’;ri3e dii-o-otion 

functions used in dovolopment 

circulation increment of vortex in uvo— dimensional flov, 
feet squared per second 

dimensionless circulation fa.ctors v i 

^v,c/<yv) 

total circulation of vortex in tv/’o— dimons ional flow, 
feet squared per second (^v,A ^v,3 + ^v,c) 

dimonsionloss la.toral coordinate of midpoint of 
specific vortex (y/s) 

circulation of specific vortex in threo-dlmcnsional 
flov, foot squared' per second 


6 


■\e 


b' 


y’ 


Bk^ Bl 


Ae, Al 


6 


F 


G 






ar 

G(cp) 


HACA TF ITo . I’l. 


swoop anglo of loading edge, positivo for swocpbocic, 
dogreos 


Symbols Pertain ing to th.o Irattorporl Method 


djjaonslonloss span along tho C.^pc lino (b/co cos A) 

dimensionless semispan along the 0.25c line (b'/C) 

dimensionless coordinate of control point along 
line parallel to 0 . 25 c line { j/oq cos A) 

dimensionless coordinate of vortex element along 
0.250c line (y/co cos A' 


perpendicular 
divided by 


distance 

Co 


from 0.25c 


line to control noinu 


distance along 0.25c line from center section to base 
of perpendicular . to control point divided by oo 


Ah - y’ 

COS“^ y'/ 3 * 



+ : ~ cos cp 

! Bh \ ^ 



/l -hi -- 1 CCS cp - C-f ) I 

V ! -Dxj o I 


■” Bp- + 


Bp a/ 


Symbols Pertaining to the We-issinger Method 
local aspect ra.tio (b/c) 

dimensionless circvilation (r(y)/bV) a continuous 
function of 
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in 

M 

n 

V 

M- 

1^1 

cpv 

'Pa 

Gv 

cv 

arv 


munber of sta.tions at ■which specific circulation 
is to be determined and at which dowmrash is 
suamed 

number of stations at which fjj .and La(v,a) 
is to be determined ’ 

denotes at which of m points specific circulabion 
ordinate occurs 

denotes a.t ijhich cf m points downwad') is summed 

denotes at which of M 'ooints n or f-p n 
I'AA^^Ai ordinate occurs 

deno'bes which of m te.'mns in interpolation function 

^'n.,A 

circular coord'inate of point n (mt/m+l = cos“'^ ri') 

circular coord:’ nate of point v (Vjt/m+1 = cos~^ r, ) 

circular coordinate cf point p (pjt/M+1 = gos"'^ •q) 

d’jnensionless. circulation at spanwise . station cpr^ 
dimensionless circuln.tion at spanwise station 


cc_ 


rn = 'Pv 
chord e.t spanwise station cpy 
specific local aspect ratio (b/cy) 

influence function 


Wh.h) 

La(v,u) 

B V n i V ^ n > -^’"v I 

^■v.,n ^ , 'V 

• > functions used in Eamhematical development 

^.vGgy.n'Sv.v^ * 

tjT(t);K j 


^’n.A’^’n,!J- 

3v 


interpolation function used in maubematical solution 
l/arv 
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DESCBIPTION OF WHIGS HTVESTIGATED 

The five wings to which the methods have been applied were 
produced by sweeping the wing panels from an existing airplane to 
five angles of sweep. Each wing then consisted of a center section, 
the two main panels, and the wo tip sections. The airfoil sections 
of the root and tip were generated by direct extension of the surface 
of the panels. The geometric characteristics of tho wings are as 
follows : 



X, 

AS ' 

(effective) 

S 

-45.2° 

0.376 

2.99 

335.5 ft 2 

-29.6° 

.405 

4.45 

282,3 ft2 

.9° 

.542 

4.47 

201.8 fts 

^1.0'^ 

.442 

4.66 

288.4 ft2 

46.4° 

i .418 

3.45 

309.5 ft2 i 


In all applications presented herein it was assumed that all 
section lift-curve slopes wore O.IO3 per degree, the average value 
of this parameter for the sections at the ends of the unswept wing 
panel. Actually, tho local section slope varied from root to tip; 
however, because of the nature of tho sections genei'atcd by 
extending the wing panels, exact values of this function could not 
be determined. Corrections to the theoretical methods to account 
for such a variation wore omitted from the computations, althou^ 
the effects of such an omission are discussed later. 

Tho loading on a wing can bo sopara.ted into tho basic loading 
(that existing at zero over— all lift) which is a function of twist, 
camber, flap deflection, and plan form; and additional loading, 
iiThich is a. function of plan form and angle of attack. For purposes 
of analysis in this report, attention has been directed solely 
toward the additional loading. The wings experimentally investigated 


2 To agree with the definition of sweep used in the theoretical 
methods of span loading prediction, sweep has been i-eferred to 
the sweep of the line joining tho quarter-chord points at root 
and tip. 
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wcrd esGcntially devoid of an,y cam'bor 
■bacic loading shown by experiiicnt tos 
ct’rvcs used a.s a basis of coiaporison. 
tbo wing was rcplccod by a. flat plate 


or twist. Any evidence of 
ronoved from tho loading 
Tl'jus, for purposes of analysis, 
and local angles of cttaol: 


bocomo synonjTncus with over— all angles of attc.ck. Tho charactoi'— 


istics of wings having camber or twist tho variation of wl;ich is 
froo from discontinuities, however, could bo dotcrained ocuall.y 
well by any of the methods simply by using the true local angle 
of -a.tta.Gl!: (as mofisurod from the angle of zero lift) at each point 
considcrod. Further discussion of this problom is given in tho 
appendixes . 


FE0CEDUS2S 

AJ.1 methods descrilbed horoin ave ext*:, ns ions of simplified 
wing, theory and so oro subject to tho same a.ssui;iptions . 


1 . Tl'ie fluid is incomprossiblo . 

2. TliO flow is potontiaJ. . 

f. Tile circulation is such that, after Kutta-TJo’alcovski, the 
stagn-ation point occurs at tho traiiling edge of the airfoil. 

4. Tlie wing is roprosontod by a. thin vortex sheet in 'ohe chord 
piano having a. plan foim idontioal to the wing plan I’crm. 

5- Aill vertical di 3 pla.nem-ont.s can bo ignoa’ed. This mo-ans, 
for insta.nce, that (a) when camber is introduced, the cherdwise 
angular variation is considoi’od but not the chordwiso vortical 
displaccmont; (b; when angle of attack is considered no vortical 
chordwiso displacements arc considered; and (c) the trailing vortex 
sheot lies always in the sane horizontal j^lnne ns tho wing. 'lijis 
ass’jmption strictly limits tho -ana-lysis to uncambored wings r.t 
zero angle of a.tta.ck; such limitations, ’lowevor, cm be mcdoratcly 
exceeded . 

In replacing the wing by a vortex sheet, t!io strength of the 
vorticity 7 at any point is related to tho difforontial pressure 
Ap at tha.t point by 


Ap = pVy 
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Tho protlom of obtaining tho loading, or distribution of Ap, 
over tho ving is thus rosolvod into that of obtaining tho strength 
of Torticity 7 witliin tho plan form. Tho control condition v’rich 
is enfoi’cod to obtain tho distribution of 7 is that no flow can 
occur through tho vorticity sheet, or in othor words, that the 
do'.mwa.sh produced by tho vorticity is proportional to tbo slopo of 
tho sheet at any point within its limits, Tiie determination of 7 
woiild bo exact if its distribution were considered continuous end 
if tho foregoing cenditien wore enforced at an infinite n'jmbor of 
points. Such an exact determination is impractical j consoquontls', 
simplifying approximations must bo introduced. Tho simplifications 
generally used are thoso cf (ij concentreting cr restricting II10 
continuous vorticity chordwiso and/or spanwise in ordei’ to make 
tho dotormination of its distr’bution amenablo to mathoinatical 
troatmont; and (2) roprosenting the distribution of vorti'dty or 
of circ'ui-ation by a mathematical expression, us^ually a sex-ios, 
containing a. finite number cf unkno^-/n coefficients whoi’e an infinito 
mrnber arc genorally required for oxact’ioss; and (3) limiting the 
number of control points at which tho condition of no flow. through 
tho siieot is sa.tisfiod. Tho differences in the various mothods 
developed for predictiag the distx’ibution of vorticity arise, 
thcroforo, fi'om (l) the manner of covicontrating or restriouing tho 
vorticity; (2) tho differences in the foian of tho mathoriiatxcal 
cxi^ressions usod to describe tho vorticity distributions; and 
(3) tho choice in number a.na. location of tho- control points a.nd tho 
precise methoxn-atical procedure used to obtain a solution. 


The F-allmcr Method ■ 


The wing ic first considered as a ontinucu.s shoot of 
vorticity whoso strength distribution is oxeprossod by the double 
scries: 


j = 8s..XjtC:A-a JT:^ 

^ -Of 

i cot ^ tao,o 

'QS-0 1 -i- ^0 ^ 2 • • • ' 


+ sin 6 (ai^o 

+ h"i,i + rjS ci^2. . .) 


+ sin 20 ( 62,0 

+ rfo^ + q ^ &r> o - • 


+ . . . STTi,n ^ 

sin m 9 

in which q = § B 

' cos -1 ^ 
c/2 



( 1 ) 
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Evaluation of 'tiiG unlcnowns am,n is porformod in tho follovr— 
ing manner by 

1. Concentrating the vorticity both chordwise end apenviso 

into a ayatem of 84 finite horaeshco vortices (fig. 1(a)) 

2 . Expressing tho circulation of these vortices in terms 

of the unicnOTOS in equation (l) (appendix A) 

3. Summing at a n'jmbor of control points on tlic wing the 

downi-ra.ah produced by all tho vortices of the subject 
system and computed by moans of the 3iot-Ssvart lav 

4 . Equating the downwash angle thus determined to the slopo 

' of the plate at those points thereby i^’orming equations 
involving tho unkno-vm coefficients 

5. Solving those equations a imulta.no ously to evaluate tho 

coefficionts 

Substitution of those values in equa.tion (l) gives the desired 
oj^ression for the load distribution. 


'Eae Ilutterperl Method 

Mutterperl considered only spanwise distribution of voi-ticity. 
In such an s.pprcacb the chordwlse distribution of vorticity is 
concentrated into the circulation of a. lifting line. (See fig. l(b) 
The distribi’tion of this circulation along the line is then repre- 
sented by the Foui’ier series . 


r = 4rtVcQ sin a 


a^n+i sin (2n-i-l) q) 


( 2 ) 


n=0 


The vinhno-vms to be evaluated to obtain tlie distribution of f 
are the coefficients o-sn+i" downwash produced at points on 

the wing by the lifting line a.nd its trailing vertex system can be 
expressed in terms of these unimoims by applica.tion of the 3iot-Sa.va. 
law to this equation. (See appendix B.) Equating the expression 
for downwash angle to tlie slopes of the mean camber lines at these 
points produces a. set of equations which contain the unhnowns a.£n+i 
simultaneous solution of these equations evaluates the coefficients. 
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The Weissinger I'fethod 

From extensions of the Multhop procedures,, Weissinser de’reloped 
two methods of obtaining span loading, one ba.3sd on lifii-ing surface 
concepts, the other on lifting line. Tne lifting Girrface method, however 
amounted to little more than a substitution of the theoretical 
additional chordwise loading, represented bj 7 = constant x Y cot 
for tile concentrated load of the lifting— line method. According 
to Neissinger the surface mothod pi-oved tc be ccnsidorably longer, 
and gave results with an accuracy only slightly superior to those of 
the line method. Fci* this reason, only the latter is described 
herein. 

As in the Mutterperl mothod, the contimious chordwise distri- 
bution of vorticity is concontratod into the circtilation of a 
lifting line. (See figure 1(c).) Tno distribution of this circula- 
tion is then specified by 

G-(cp) = ^ Gn ^ 3in PiCpn sin Picp (3) 

n=l (^1=1 

Tho circulation F(y) is represented nondimensionally as G(p) 
in this expression and the unknov/na to be evaluated are G-^, the 
circulations at specified locations along the line. The dot-mwash 
produced at points within the plan form by the lifting line and 
its trailing vortex systorn can bo expressed in terms of these 
unknowns through application of the Biot— Savart law to equation (3). 

(See appendix C.) liquating the expressions for doimwash angle so 
obtained to the slopes of the mean cambex- lines at these points 
results in a set of equations with oininovms Gn. Simultaneous 
solution of these equations evaluates tho unkno\ms. 


EXPSP.LMEKTiiL DAT.\ 


Pressiire dam vrere taken at a tunnel speed of 9G mi3.es per 
hour 'vaich corresponds to a P.eynolds number cf approximately 
9,000,000. Data were talcen over an angle— of— attack range —3° to 
9®. Plots of the chordwise distribution of prossiure coefficient 
P =(p— ?Q)/q at several spanwise stations were dra’./n and integrated 
to obtain the local lift at these stations. Tliese values of local 
lift were then plotted against angle of attack, and tlie resulting 
local lift— curve slopes were used to obta,in the curves of the spanwise 
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distribution of additional load and of additional lift coefficient 
shown herein. The maximum ori-or in any loca.1 lift— curve slope as 
the result of scatter^ etc., is estimp.ted to he 0.C02 per degree. 
Such an error would produce a variation of the distribution curves 
of about one-ha.lf to one— third the magnitude of the discrepancy 
between the theoretically computed and the oxperimon bally obtained 
curves . 


KSSULTS AI© riSCUSSION 

Comparable spanvrise distributions of the loading coefficient 
cZCi^/OLCav as calculated by the three methods and as detenained 
from the oxpc-rimental surveys are presented in figure 2. Similar 
presentations of local lift coefficient cj/Cp are presented in 
figure 3« theoretically predicted values of lift— curve slope 

dCL/do, and spanwi.^se centei- of pressure yjosition for the different 
wings ara presented in table I. 

From figures 2 and 3 it is apparent th<at all the theoroticai. 
methods tend to predict higher loadings at the center and lower 
i.oadings at the tip than were moasured. In general, the Fallmer 
method errs less in this respect than do the others. Mutterpierl 
distribution repireoontationo for the swept— back and xmswept wings 
are only slightly less accurate than those of Fslkner. On the 
other hand, for the swept— for\/ard wings the Muttcrporl distribu- 
tions departed from the experimental distributions to the extent 
that they must bo considered unusable. Weissinger distribution 
representations wero equally accurate for svrept— back and swept— 
forward wings. The average accuracy for this method was only 
sligbtlj^ less than that of the Falknor method. 

In regard to the center of pressure position and lift— curve 
slope, the closest predictions in all instances were those made by 
the F-a.lknor method. The Muttorpcrl method, in the range in which 
its applications may be considered usable, was also quite QCCUi’o.to. 
The Weissinger method gave good center— of— pressure positions in all 
instances and accurate values of lift— curve slopes in all instances 
except for the swept wing. 


The time studies of the calculations indicate that the Falkner 
method takes from 2k to 37. hours. The groator part of this time, 

16 to 20 hours, is consumed in determining the values of the down- 
wash factor F for the different vortices. Uie major part of the 
remaindor is needed for the solution of the simultaneous equations, 
which often prove to bo ill conditionod. Tho Muttcrporl method 
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takes from. 20 to 28 hours, the greater part of the time oeing 
consumed in the Siiopaon mle integration of the factors Fg’ to Fg* . 
The VTeiss.inger method using m = M = "i , takes only 2^ to 3 hcvirs, 
in which there is no phase that consumes an outstanding amount of 
time . 

It has been stated prev.iously that the section lift-~curve 
slope c-i^ of all sections on all five wings was assumed to ho 
0.1030 per de£p.'’ee . Tiie thiclmess variations from root to tip, 
however, iiidicate that variations in prohahl.v exist for each 

wing. Unfortunately, tlie distortions of the sections resu.lting 
from the manner in which the wings vrere constructed preclude an 
exact determination of what the variation migiit he for all hut the 
unswept wing . For this z’eason, the readily applicahle corx’ection 
to theory (see appendixes) for a variation in ci^ \ms not 
included in the computations. TOiile this correction would 
account at least in part for the aforementioned discrepancies 
between theoretical and experimental loading distributions, it 
should net alter the relative evaluation of the three methods. 

In considering the three methods it should he noted that two 
of them, those of Weissinger and Mutterperl, have identical aero- 
dynamic a.pproaches and differ only in the mathematical procedure. 

It would he expected, therefore, if no compromise were made in the 
mathematical accuracy (i.e., if a. large n'omher of terms were used 
in the series), identical results would he obtained. Further, if 
similar limita.tions were impressed upon the two methods it might 
well he a.ssimied that results of comparable accuracy would ha 
obtained. The failux'e of the Mutte.rpei'l method to predict 
acceptable loadings on the swept— f orward wings is inexplicable on 
these grounds and, as a result, must he attributed to an inconsist- 
ency introduced in the mathematical development. An additional 
advantage of the Weissinger method is that it lends itself to the 
pre tabulation of a number of constants which ai'-e applicable to Phe 
solution for any plan form. It is because of this that the Weissinger 
metliod proved loss tiaie consuming than that of Mutterperl which 
cannot be handled in this manner. In general, then, it is apparent 
that the Weissinger method offers sevez'al advantages ovei- the 
Mutterperl method, which, however, stom entirely from the mathematical 
pihase of the solution. Insofar as the aerodyicamic concepts are 
ooncomed neither method should be expected to be superior. 
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Tlie Nal-kncr method offers a definite aerod'-naraic a.dvantago in 
bbat the wing is represented hy a lifting surface r<athrr than a 
lifting line. Fa-om a oonsidaration cf only the spanwise distribu- 
tion of loading, the time required to use the Fallmci’ method appears 
excessive when the very mmior impro’rement in accuracy is rccognizod. 
However, if surface load.ing or chordwise loading wore desired, the 
method would imdoubtcdly show marked superiority. The relatively 
long period, of t:Lme required to obtain a solution by this method 
is in great mca.3iu'’e 0. result of the largo number of purely mcclir.nlcal 
functions inherent in the- method. It can be cjcoectcd that suc.h 
processes are amenable to handling by mechanical means if smTicient 
use is to be made of t.hc method to won-ant their construction. One 
such aid of relatively simple foi-m has been a.rjplied in other sp"-n~ 
loadi.ng computations using the Fal’mcr moukod and i-csulted in cutting 
the computing time by 30 percent with no serious loss in accuracy. 

It consisted of constructing a. large-scale contou.r ch.ai’t of the do";vn— 
wash field ai-ound, a horseshoe vortex and using this in conjunction 
with an appropriately sealed dj-awing of the wing to rood directly 
the dcvmwaoh at the various control points. 

A further adveantago of the Fallcncr method over the lifting—l'hcc 
methods can be seen in the' increased flexibility resulting from the 
system of finite vortices which permit application of tlr^s method uo 
a variety of plan forms beyond the scope cf the other methods the 
lifting lino pattern and oontrcl—point positions of whies arc 
fairly rigidly specified. In this re.^i.rd, Ihalloacr 2'ias successfully 
aqjplicd tho method to a. pterodactyl wing and to a wing with a. 
parabolic 0.2pc line. It should bo i-onombercd, hovmvcr, that 
should the plan, fcrri be of su':.h a nature .as to require a Ecdifica— 
ticn cf the vortex lattice, tl:c work iiivolvcd will be considerably 
increased. 


COIICLtJSIOIiS 

From the results of the subjece investigation the foilowin.g 
conclusions have been drawn: 

1. Wixere an over— all study of the effects of sweep and pla.n 
form on span leading, lift— curve slope, etc., is desired a.nd 
where good accuracy is desired for minimu;.i e.rfcrt, tho Noissingor 
method is most useful. 


16 


MCA 'M Wo. 


2. Where a detailed study of a. specific wing ?s desired and 
utmost o.ccuracy is important oven at the expense of considerehle 
computing effort^ the Falkner method should be used. 

3. The Muttcrperl method offered no a.dvantagcs over the o'chcr 
methods either in teians of accuracy or facility. 


Ames Aeronautical Laberatory, 

Watlonel Advisory Committee for Aeronautics, 
Moffett Field, Calif . 
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APPSITTIX A.- PEETII'ieilT’ INFOm^TIOK AMD COME’UTIEG PEOCEDUEES 
FOE USE WITH THE F.UEG7E:E METHOD 


Selection of the Vortex Pattern 


The relative strength of the circulation of the vortices in a. 
network, as expressed in terms of the unlcnowns in the series 
equation (i), depends upon the vortex pattern and the terms in the 
series only, not upon v/Ing shape. Tables of such circulations 
can be set up for use with any specified pattern. Faikner, on tlio 
basis of ills applications, selected the pattem of 84 vortices 
shown in figui’e 1(a) as suitable for raost wings. Wnile it is 
recognized that other patterns mi-ght produce more accurate results 
In particular instances, the advantages of uhis regular patteiTi 
in reducing the computational work are gr*eat and hence it was 
used for all applications included herein. 


Limitations of the Series 


Tlie number of terms in equation (l) required to obtain a good 
approximation of the load distribution depends on the rapidity with 
which the series converges for each a.ppli.catio;o. For the calculation 
of symmetrical loading, Faikner concluded that a minimum of t}:ree 
chordwise and three spanwise terms (nine unkno'.inr'.s) should be used 
for all swept wings, while a minimum uf two sparwise and three chord- 
wise terms (six unlocovms) should be used for straight wings. 

It should be recognised tha.t, a.s it is given, this aeries 
will not converge when attempting to approximate a sui’face leading 
where discontinuities exTst such as those resulting from partial 
span flaps. A slii^it modification to the series, however, will 
enable it to approximate the loading v/hore such a discontinuity 
occurs. Faikner has determined tne necessary modification in Ills 
investigation of wings with flaps and ailerons deflected. 


Determination of Circulation of IJotwcrk Vortices 

Ones the vortex pattern and number of terms in the basic 
series have been established, the circulation of the vortices a.s 
expressed in terms of the unioiowrs in equation (l) can be determined 
by replacing the continuous vorticity chordwTse and. spanwise of 
equation (l) with the concentrated stepped leading of the lattice. 
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The chordviGO concendraticn of the load is detcnii'j.ned the 
conditjon that^ ab points located Eild.¥ay between the loads (ao 
one-quarter, one-half, ard three-quarte'-- chord), the dovnwash^ 
produced by the four chordwise load.s Tv be the same as wouxd be 
produced by the continuous chordwise loading in two-dimensional 
flow, and the limitation that the sum of the isolated loads be 
equal to the integral of the continuous load. 

When only the first term in the chordwise series of equation 
( 1 ) is considered, 


7 


8sV ba n a 
c 


00 



n=0 


(Al) 


or, since only chordwise loading is being considered 
which are not a function of the chordwise variable 


e 


into a constant A where 


all factors 
can bo put 


x-nen 


= 8s tan a ij l-q‘"'^ 


r i_ T 

n=0 




(A2) 




(A5) 


Then it can be shown that if the flow is considered two dimens_onc:.l 
tho . downwash angle at any point along the chord is 


w _ 1 A 
V 2 c 


/ ydc = W j cot |- sin 9 d0 

J-0/2 


flAV 

o 


(a4) 


and 
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Tlio following foiir equations may therefore he obtained: 





whoro equa.tion (a) equa/bes the into grad of the continuous loading 
to the simi of uho circulations of the separate loads^ and oqiiations 
(h), (c), and (d) equate the dow.easb au the throe pivotal points 
as produced by the continuous loading to that produced by tho fom’ 
loads of Fv. A simultaneous aol’ation of those equations gives the 
increments of circulation of the four cherdviso vortices which are 
equivalent to the continuous loading roprosentod by terra 1 of 
equation (l) . 

A simiJ.ar solution when 



wnorc 


00 



n=0 
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••• 


(A8) 


or if 



= AV;tGv ^ 


+ BVrtGy g + CVTTGa; q . . . 


(A 9 ) 


The suhvstitution into this equation of the values of k, '£>, etc.^ 
introduces the spanwise variable r, . 


Since the circiilation of specific horseshoe vortices is now 
being considered, the circulation no longer varies continuously 
■along the span but reMains constant tiiroughout the length of the 
bounded lines. Ibis is equivalent to the assumption that the 
continuous loading is stepped at intervals equal to the length 
of the bounded lines of the network vortices. The continuous 
variable tT of equation (l) or equation (AlO) is therefore replaced 
in a new equation by specific values p of which indicate the 
midpoints of these lines. This new equation which expresses the 
circul-ation of any network vortex is then 


"'The values of Gv^a Gp^p and Gy ,C presented by Falkner in 
reference 1 were found to be in error. Under the direction 
of Mr. Arthur Jones these valiies were recomp”.ted at Areas, 
and the values so obtained are presented in table Al. 


rv,n = StcsV tan a [Gy, A i^o,o + h ^c,x + ao,2---) 

+ Gv,e (ai,o + p e.1,1 + • •) 

+ Gv,C (a-2,0 + 'iT ©2,1 + + •••) 


+ . . . ] 


(AlO) 
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Tv.u = tan a/l-^2 + m + ...) 

+ S,B ^ + ...) 

+ S,C (^2,0 •*• ^2,1 •" ^2,2 + •••) 

+ ••' ^ (All) 

or for a symmetrically loaded wing 

~ 8rtsV tan a^^l— la*- [Gy^ys '.ao^o '*" ^Oj 4 ■*■ •••) 

+ S,B K,0 •*■ ^1,2 + ®1,4 + •••) 

+ ^VjO (^ 2^0 + 3 - 2^2 + ^ 2_,4 + •••) 

+ • ”^ (All)’ 

Examination of this eqioation will show that, as has prevTcisly hasn 
indicated, the l-mo-.m parts of the equation p, Gv,a/ G v,B, etc., 
are indep endent of wing shape. The products of these factors 
p v'l-l-^^ Gv,A^ d^ v'^1 — G-v,3j ■ etc., have been tabi/lated for -’vSe in 
anj”- application in which the 34 — vortex pe.toern is used. (See table All. 


Selection of Control Points 

Since one equation is fcrned at each point and since there 
should be the same n’mnber of equations as there are unknowns, the 
total number of points selected is determined by the total number of 
unlaio-wns reta,ined in the series equation. Further, the spanwise yuri 
chordwise distribution of control polntc must correspond to the 
number of spanwise and chordwise terms retained in the series. Tlie 
locations of the points chordwise and spanwise are limited to nosi— 
tions midway between or on the center line of the vortices. .Aside 
from these limitations, the exact choice of location remains a matter 
of experience. Falkner found that for a calculation of symmetrical 
loading the a,rrangoment presented in figure 1(a) resulted in good 
accvi’acy for wings with sweep. This arrangement has been used in 
all the calculations presented herein. 
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tion r 


Determination of the Downwa.sii 

downwaoh produced ty a. simple horseshoe vertex of circula— 
is expressed hy (refei-’ence 1) 


w = r 

V hr. yv V 


(A13) 


e'here 


x^'- -v(x-^)2 + (y--H-l)2 
x-'^(y^^+l) 


The downwash produced by a network vortex is then from equations (All), 
(A13), and (A14), 


x-'~y(x^)2 + (y»-l)g 

x^n'y^^l) 


(Al4) 


^ = 40 tan a •/ l—u ^ [Gv^a (^OjO + •••) 

GVjB (s-i.o + '4 si.i + 

+ Gv^C (®- 2 ,o + a2,i + •••) 

+ ... ] F (A15) 


or for a siTimictrical wing 

^ = kO tan aVl-ti^( Gv,a (a-o,o + a-o,2 + •••) 

+ Gv,B (ai_,o + ai ^2 + ...) 

= Gv^C (S2,0 + £>2,2 + •••) 

+ ... ] F (Al6) 

The coordinates x*, j""*, and consequently the factor F con he 
determined readily from wing geometry. In this regard, plots of the 
function F versus x* from 0 to 20 have been prepared at values 
of y* =0, 2, k, 6 ... 40j however because of tneir size these 
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charts as such havs not "been included in the report j but the tabular 
da.ta necessarj'' for their construction are ,gi.ven in table AIII. In 
addition, examination of (Alii) will reveal that if y* is constant 


For this reason table AIII contains onl/ positive values of x->^, 
and the function Fi + F2 is presented in table AlY. 

Summation of the domwash at any control point now results in 
an expression containing the unknowns Oj^ ^ and their numeri cal 
coefficients vrhich are products of the tati’ulated values 
G-v A.J regard it should be noted that in 

thv 3 suMmd cions'’ for a s;unmotrical (about root chord) wing, the down- 
wash factors F for symmetric ally located vortices may be added 
together prior to tho multiplication of these factors by the circula- 
tions of the vortices, sixice in this instance the circtilation of suah 
a pair of vortices will be identical. 


To obtain the additional loading, the wi.ng is considered .a flat 
plate the slope at which any point is tan a = y. Substitution of 
chis value into the do-imwash expi'easions, as evaluated at the 
several control poiirts, results in a set of equations vioh 'onlarowns 
am n* Simultaneous solution of those eq'uations eval'uates the 
'uhkno'wns which can then be introduced in'co equation (l) to produce 
the desired expression for additional loading. 

The following expressions can now bo derived readily from 
equation ( 1 ) : 


F (-x ->0 = Fi + F2 - F (x*) 


(AIT) 


whei'e Fi + F2 is a function of only 



(AlS) 


Solution for Additional Loading 



dCp _ AE 
da 16 


da 16 
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£L.£IL- = 32yi-T,^ [2ao.o4ai.o+ ^Y^(2so.g+c, 

Cjj Cp.-V ^(l6ao^O+69i,0+^2!.o^2^^^1^2"^^^0^4'*'^1^4) 

7, = 32 [3$(Sao.o+ri.c)^+l4(2ao.g+ax.g)+8(25Q^4+ax^4 )]^ (A22‘) 

^ 10^it(l6aQ^ Q+8a.j^^ Q+4fv ^2+ 2a.j^^2+2u-Q^_j+a.j^^_^ ) 

into which the coefficients Cjjj must be substituted to obtain 
the quantities indicated. 


Solution for Basic Loading 

The determination of the basic loa.ding on a wing with camber 
and twist can be ancoraplislied in seo-eral i-niys. The simplest of these 
is to calculate the total lcadj.ng, basic plus additional, at some 
finite lift coefficient and then to subtract from this the additional 
loading as calculated for tl;at lift coefficient. A solution for the 
total loading on c. ca.mboi''od and twisted vring is identical with that 
of a flat— plate wing up to the formation of the 3im'ul.tc;neou3 equations. 
For the flat— plate wing all local geometric angles of attack were 
identical to the wing geometric angle of attack; in this instance 
local goomoti’ical angles of autack are in addition a function of 
the cambor and twist. 

If the midwing section of the wing is chosen o.s a reference and 
sot arbitrarily at some angle ag, then bbo local geomotricai. 
angles of attack at the various control points are l-cno''-m exactly; 
however, tlae angle a of the reference from the zero ].ift angle of 
the wing is not know:a. To obtain a solution under those circumstances 
the values tan aiocai and tan as are substituted for w/V and 
tana, respectively, in the downvra.sh exioressions, and a solution for 
the coefficients a;m,n obtained in which, however, these coeffi- 
cients will ho in error by the factor tan a/tan ag. If those coeffi- 
cients and the factor tan a are then introduced into the coqorossion 
for lift coefficient. 


L = tan a ( l6a 0 ^ 1 ^ 0+^- c , a+^a. ^ ^ g+Sa.^ ^ ^ J (. .23 ) 


tho lift coofficiont for the angle of attack a vill bu obtained^ 
since the error introduced by using cxg -will bo negated by tho error 
in tho coefficients Cm^n- other words, the result is the same 
as if tho correct values of ain^n boon inserted into 
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eqmtion (A23). Similarly introduction of tan., ag and the incorrect 
values o.f into the following: 

CjCt, = 4tto tana*/l-T)2 (28.^^^+a^^J] 

(A2h) 

will result in the values of the ordinates cf the total Inading 
curve for tan a Now if a solution is effected for the additional 
loading^ as previously descrihed^ and the value of the lift— curve slope 
dCL/da tlius obtained from equation (A19) is divided into the valde-'' 
of Cx, obtained from equation (A 23 ')j the correct value of the wing 
angle of attack tan a will result. If this value and the coeffi- 
cients 3 ^ of the additional loading are then substituted into 
expression (A24)^ the ordinates of the a,dditional loading curve will 
be obtained. Subtraction of these from the ordinates of the total 
loading curve will result in the ordinates of the desired basic 
loading curve. 


Correction for Section Lift— Curve Slope 

Through the general development of the method all section lif1>- 
cuerve slopes were assumed to be the theoretical 2n per radians 
(0.1096 per deg). As this assumption is not valid for aJ.1 sections 
the final expression for vorticity will be in error unless a correc- 
tion is applied. If the section lift-curve elope is the same at all 
sections of the wing, the error may be corrected by simply multiply- 
ing each coefficient ^ ratio of actual section lift- 

curve slope to tliaoretical section lift-curve slope. A varying 
section lift-curve slope can be accommodated almost as easil 5 '; 
however, in this instance the correction must be introduced into 
equation (l) as a function of the spanwise variable ?[. 


Computing Instructions 

The following instructions apply to unyawed straight tapered 
svrspt wings without camber or twist. 

The coordinates x*’ and y*' relating all vortices and 
control points to the center section leading edge of the wing are 
calculated on form A(l) using the relations 
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= 20 |i 

X*' =20 I ^ ! tan A - [1- (l-X.) |n|] 

[jL-rX }RR 


for the vortices^ and 


(A25) 


X*’ 


= 20 T] 

= 20 in I tsnA - 

t : I 


80v 


[1 - (1-X) jn|] 


(A26) 


for the control points. Since the wing is syinaetrical x*' will he 
the same for siniilarlv' loce.ted vorticos on each wing half, and the 
values of y-' for the left wing will he the same as those for the 
riglit wing, although of opposite sign. For this reason these values 
only need he computed for positive values of or p. 


The values of x^*- and relating a control point to each of 

the vortices are obtained hy subtracting the values of X'"’ and 
y^-’ of the vorticos from those of the control point, column 9 3^ 

form A(1) from a value in column l3 or 12, form A(l), respectively; 
X* and y’'" are then tahuLated on a form A(2), using a separate 
form for each control point. It should he noted that since the 
coordinates of the vortices at p = O .9625 are based upon a unit 
length y^ one— quarter normal size, x'*^' and y* for these 
vortices are four times the normally calculated values. Lastly, by 
virtue of symmetry of plan form, the coordinates can he tabulated 
so tha.t two values and yp>* exist for every value of x’*’. 

These coordinates are now used to enter charts of the downw8,sh 
function F ••as prepared fi'om the values in table AIII. The values 
obtained for the vortices at p = O .9625 should be multiplied by 
four. Booauso of symmetry of loading, Fj^ and Fp_ can be and are 
added together. 

Tlie simultaneous equations set up in tabular form in form A(3; 
are now obta.ined as follows: Considering the first equation or 

column 1, the second number, the numerical coefficient of ao.Oj 
is obtained by multiplying the values of Fp + Fj_ in coluimi t, 
form A(2) as determined for control point 1 by the values listed 
under Uq q in table AH end simnming the products. Similarly, 
the third ■^number in col’Emi 1, the mmiorical coefficient of •3.x,o-’ 
is obtained by m>oltiplying the values in col’usmi 7 by the values 
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listed under in tatle All and suinming the product. The process 

is repeated using the values listed under etc,, until 

the entire equation is ohta.ir.ed. 

The second equation, column 2, form A(3), is set up in the same 
manner except that the values in form A(2) as determined for a 
second control point are used. The procedure is then repeated 
until the nine i-oquirod equations arc formed. 

The constant numbers, row 1 of form A(3)^ Qfc obtained as 
follows: The dewnwash at control point 1 is 


In like manner the constants for the other equations are also 

0.0250, 

Tlie equations are set up in this manner to facilitate their 
solution by the method outlined in reference 4. Of the various 
methods for solving a large number of simultanc-ous equations 
which were tried, the method of reference 4 w.as found to be most 
rapid and st.raightf orward where only standard computing ma.chines 
were a.vaila.ble , 


— = 40 tan a X column 1, form A(S) 


or 



Equating ~ to the slope of the p.late, tan a 


40 


= 0,0250 = column 1, form A(3) 
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Form f \-~ Computing form for Falkners method 

(^Underscored numbers are sample calculations) 


( 3 .) 
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TABLE AI.- CIlORDVfISE FACTORS G FOR VORTEX 
PATTERN UTILIZING FOUR CHCRDV«SE VORTICES 


V 

<^V,A 

*^V,B 

^v,c 

0.125 

0.27337 

0.04902 

0.07282 

.375 

.11630 

.07593 

.03323 

.625 

.06947 

.07598 

-.03823 

.875 

.04056 

.04 902 

- . 072 82 
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TABLE All.- PRODUCT OP CHORDWISE AND SPACTIISE FACTORS P® am.n TO BE 
USED WITH THE 84-V®TEX PATTERN SIMJBTRICAL LQADIHO 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


V 


* 0 V,A 

for 

^ 0,0 

X Gv,B 
for 

*t ,0 

X Gv,C 
for 

^ 2,0 

X Gv,A 
for 
*^ 0,2 

X Gv.b 

tor 

X Gv.C 
for 

‘0.4 

Ov.B 

for 

ai,4 

X Oy,C 
for 

a.,4 

0*125 

0 

0.27337 

0.04902 

0.07282 

0 

0 

0 

0 

0 

0 


• 1 

.27200 

.04877 

.07246 

.00271 

.00049 

.00072 

.00003 

0 

.00601 


•2 

,26786 

.04803 

.0715$ 

.01072 

.00192 

•00285 

.00044 

.00005 

.00012 


•3 

.26077 

.04676 

.06947 

.02348 

.00421 

.00626 

.00210 

.00036 

•00066 


•4 

.26054 

• 04493 

.06674 

.04008 

.00719 

•01068 

.00642 

.00115 

.00171 


o 6 

.23674 

.04246 

•06306 

.06918 

.01061 

.01677 

.01479 

.00266 

• 00394 


• 6 

.21870 

.03922 

.05826 

.07873 

•01412 

.02097 

.02836 

.00608 

•00766 


•7 

.19621 

.03601 

.05200 

.09666 

.01716 

.02648 

.04688 

.00841 

.01249 


.8 

.16402 

.02941 

.04369 

.10497 

.01882 

•02796 

.06719 

.01206 

.01790 


• 9 

.11916 

.02137 

.03174 

.09663 

.01731 

.02671 

.07818 

•01402 

.02083 


.9626 

.07417 

.01330 

.01976 

.06870 

•01232 

.01830 

.06364 

.01141 

.01696 

o576 

0 

•11680 

.07698 

•03825 

0 

0 

0 

0 

0 

0 


.1 

•11622 

.07560 

•03804 

.00116 

.00076 

.00038 

.00001 

.00001 

0 


•2 

.11444 

.07446 

.03746 

.00468 

.00298 

.00160 

.00019 

.00012 

•00006 


• 3 

.11142 

.07248 

.03647 

.01003 

.00663 

.00328 

.00090 

•00059 

.00029 


•4 

.10706 

•06964 

.03604 

.01712 

.01114 

.00660 

.00274 

.00179 

.00090 


•6 

.10116 

.06680 

.03311 

.02529 

•01645 

.00828 

.00632 

.00411 

.00207 


• 6 

•09344 

.06078 

.03068 

.03364 

•02188 

.01101 

•01211 

.00788 

•00396 


.7 

.08341 

.05426 

.02730 

.04087 

.02669 

.01338 

.02003 

.01303 

.00666 


o 8 

.07008 

.04669 

.02294 

.04486 

.02918 

.01468 

.02871 

.01868 

.00940 


• 9 

.05091 

.03312 

.01666 

.04124 

•02683 

•01360 

.03340 

.02173 

.01093 


• 9626 

.03169 

.02061 

.01037 

.02 934 

.01909 

.00961 

.02709 

.01769 

•00890 

o625 

0 

.06947 

.07698 

-.03823 

0 

0 

0 

0 

0 

0 


.1 

.06912 

.07660 

-.03804 

•00069 

.00076 

-.00038 

.00001 

•00001 

0 


o 2 

.06807 

.07446 

-.03746 

.00272 

.00298 

-.00160 

•00011 

.00012 

-.00006 


• 3 

.06627 

.07248 

-.03647 

•00697 

.00663 

-.00328 

.00063 

.00069 

-.00029 


o4 

.06367 

.06964 

-.03604 

.01018 

.01114 

-.00660 

.00163 

.00179 

-.00090 


•6 

•06016 

.06680 

-.03311 

.01604 

.01645 

-.00828 

.00376 

.00411 

-.00207 


• 6 

.06668 

.06078 

-.03068 

.02001 

.02188 

-.01101 

.00720 

.00788 

-.00396 


o7 

.04961 

.05426 

-.02730 

.02431 

.02659 

-.01338 

.01191 

.01303 

-.00666 


• 8 

.04168 

.04669 

-.02 2 94 

.02668 

.02918 

-.01468 

.01708 

.01868 

-.00940 


• 9 

.03028 

•03312 

-.01666 

•02453 

•02683 

-.01360 

.01987 

.02173 

-.01093 


• 9626 

.01886 

.02061 

-.01037 

.01746 

.01909 

-.00961 

.01617 

.01769 

-.00890 

• 876 

0 

.04036 

.04902 

-.07282 

0 

0 

0 

0 

0 

0 


a 

.04016 

.04877 

-.07246 

.00040 

.00044 

-.00072 

0 

0 

-.00001 


•2 

.03954 

.04803 

-.07136 

.00168 

.00192 

-.00286 

.00006 

.00008 

-.00012 

] 

o3 

.03860 

.04676 

-.06947 

.00347 

.00421 

-.00626 

.00031 

.00038 

-.00066 

1 

o4 

.03699 

.04493 

-.06674 

.00592 

.00719 

-.01068 

.00096 

.00116 

-.00171 

I 

• 5 

.03496 

.04245 

-.06306 

.00874 

.01061 

-.01677 

.00218 

.00266 

-.00394 


• 6 

.03229 

• 03922 

-.06826 

.01162 

.01412 

-.02097 

.00419 

•00608 

-.00766 


•7 

.02 882 

.03501 

-.05200 

.01412 

•01716 

-.02548 

.00692 

.00841 

-.01249 


.8 

.02422 

.02941 

-.04369 

.01650 

.01882 

-.02796 

.00992 

.01206 

-.01790 


o9 

.01769 

.02137 

-.03174 

.01426 

.01731 

-.02671 

.01164 

.01402 

-.02083 


• 9626 

.01096 

.01330 

-.01976 

.01014 

.01232 

-.01830 

.00940 

.01141 

-.01696 1 



TABLE AIII.- DOWIWASH FACTOR F IN THE FIELD OF A HORSESHOE VORTEX 

[At positive values of x* only] 



"{ 

F + ( 2 /x) 

F 

X*' 

y-”- = 0 

CM 

li 

y* = 4 

6 

8 

10 

12 j 

11. 

i 16 

18 

20 

0 

2.00000 

—.66667 

-.13333 

-.05714 

-.03175 

-.02020 

-.01399 

-.01025 

1 

-.00785 

i 

-.00619 

-.00502 

.1 

1.90020 

-.62234 

-.12976 

-.05616 

-.03135 

-.01999 

-.01387 1 

! -.01018 

-.00780 

-.00616 

-.00499 

.2 

1.80200 

-.57875 

-.12625 

-.05519 

-.03093 

-.01980 

-.01375 1 

i -.01011 

-.00774 1 

-.00612 

-.00497 

.3 

1.70647 

-.53652 

-.12271 

-.05421 

-.03053 

-.01959 

-.01364 ! 

! -.01003 

-.00770 1 

I-. 00609 

-.00493 

.4 

1.61484 

-.49621 

-.11920 

-.05326 

-.03011 

-.01938 

-.01352 1 

-.00996 

-.00765 

1-. 00605 

-.00491 

.5 

1.92788 

-.45819 

-.11573 

-.05226 

-.02973 

-.01919 

-.01340 1 

1 -.00988 

-.00760 

-.00602 

-.00489 

.6 

1 .44603 

-.42269 

-.11228 

-.05130 

-.02934 

-.01897 

-.01329 1 

-.00981 

-.00755 

-.00599 

-.00486 

.7 

1.36954 

-.38981 

-.10889 

-.05034 

-.02893 

-.01878 

-.01317 

i -.00974 

-.00750 

-.00595 

-.00484 

.8 

1.29844 

-.35957 

-.10555 

-.04938 

-.02854 

-.01857 

-.01305 

j -.00967 

-.00745 

[-.00592 

—.00481 

.9 

1.23254 

-.33186 

-.10227 

-.04843 

-.02814 

-.01838 

-.01293 

-.00959 

-.00740 i 

1 -. 00589 

-.00478 

1.0 

1.17158 

—.30655 

-.09904 

-.04750 

-.02774 

-.01817 

-.01282 

j -.00952 

[-.00736 

-.00584 

-.00476 

1.5 

.92963 

-.21017 

-.08400 

-.04292 

-.02581 

-.01718 

-.01224 j 

i -.00915 I 

-.00711 

— . 00568 

-.00464 

2.0 

.76393 

-.14956 

-.07093 

-.03863 

-.02394 

-.01620 

-.01167 

-.00880 

1 — .00686 

-.00551 

-.00451 

2.5 

.64594 

-.11032 

-.05986 

-.03468 

-.02214 

-.01526 

-.01112 

-.00844 i 

; - .00663 

-.00534 

-.00439 

3.0 

.55849 

-.08398 

— .05066 

-.03107 

-.02045 

-.01435 

-.01057 

-.00810 

-.00639 

-.00517 i 

[-.00427 

4 

.43845 

-.05255 

-.03683 

-.02492 

-.01739 

-.01264 

-.00953 

-.00743 

1 -.00593 

-.00485 

-.00403 

6 

.30574 

-.02555 

-.02100 

-.01632 

-.01255 

-.00974 

— . 00770 

-.00620 

i -.00508 i 

i-. 00423 i 

I-.00357 

8 

.23444 

-.01489 

-.01318 

-.01112 

-.00918 

-.00752 

-.00620 

-.00515 

-.00433 1 

!-. 00367 

-.00314 

10 

.19012 

-.009691 

j -.00893 

-.00791 

—.00685 

-.00587 

-.00500 

-.00428 

! -.00367 i 

i-. 00318 

-■.00277 

15 

.12889 

-.004381 

1 -.00422 

-.00397 

-.00368 

-.00335 

-.00305 

-.00274 

j -.00247 ^ 

-.00223 i 

1-. 00200 

20 

.09750 ! 

-. 00248 i 

-.00243 

-.00234 

-.00223 

-.00211 ! 

-.00199 

-.00185 j 

—.00171 1 

!-. 00157 i 

i -.00148 
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TABLE Ccncluded 


X * 

1 ' * ■ ■ 

F - • ■ V . • ‘ . 

. • ‘ i 

1 

y * 22 

24 

26 

28 

30 

, 1 

r . 

v '! i 

r . • 

; ,36 . i 

1 . i 

38 •; 1 

‘ ^ K 

0 

-.00413 

-.00348 

-.002 96 

-.00256 

-.00222 

-.00196 i 

-.00173 j 

1 ^ 
1 -.00154 j 

[ 

-.0015^ 

•-.00125 

.1 

-.00412 

-.00346 

-.002 95 

-;0C255 

-.00221 

-.C01S5 

-.00172 

; -.00154 j 

! -.00139 

; -.00125 

.2 

-.00410 

-.0034 5 

-.002 94 

-.00254 

-.00220 

-.00195 

-.00172 I 

! -.001.54 

-.00138 

1-;. 00124 

.3 

-.00409 

-.00343 

- . 002 93 

-.00253 

-.00219 

-.00194 

-.00171 i 

j -.00153 

-.C0158 

. -.00124 

.4 

-.00407 

-.00342 

-.002 SI 

-. 00252 

-.002 1 8 

-.00193 

-.00171 

! -.00153 

-.00137 

-.C0123. 

.5 

-.00405 

-.00340 

-.00290 

-. 00251 

-.00218 

-.0019-2 

-.00170 

i -.00152 

-.00137 

•-.00123 

.6 

-.00403 

-.00338 

-.00289 

-.002 5 0 

-.00217 

-.00192 

-.00170 

-.00152 

-.00137 

00125 

.7 

-.00401 

-.00337 

-.00288 

-.00249 

-.00217 

-.00191 

-.00170 

-.00151 

-.00136 

:-. 00122 

.8 

-.00400 

-.00336 

-.00287 

-.00248 

-.00216 

-.00190 

-.00169 

-.00151 

-.00136 

■-!. 00122 

.9 

-.00398 

-.00334 

-.00286 

-.00247 

-.00215 

-.00190 

-.001 88 

-.00150 

-.00136 

-.00122 

1.0 

-.00396 

-.00333 

-.00285 

-.00246 

— . Cu2 1 5 

-.00190 

-.00168 

-.00150 

-.00135 

-. 00122 .; 

1.5 

-.00387 

-.00325 

-.002 7 9 

-.00242 

-.00211 

-,.00186 

-.00166 

-.00148 

-.00134 

-.00120 ■ 

2.0 

-.00377 

-.00318 

-.002 7 3 

-.00237 

-.00208 

-.00183 

-.00163 

-.00146 

-..00132 

i-.OOLlS ■ 

2.5 

-.00368 

-.00311 

-.00268 

-,00233 

• -.00204 

-.00180 

-.00161 

.-.0Q144 

, -.00130 

-.00117 

3.0 

-.00358 

-.00304 

-.00263 

-.002 2 8 

-.00200 

-.,0017.7 

00158 

-.05842 

■ -.00128 

-.00115 

4 

-.00339 

-.00291 

-.00251 

-.002,19 

-.00193 

-.00172 

-.00153 

-.00137 

-.00124 

-.00113 

6 

-.00305 

-. 00263 

-.00230 

-. 002 ; 0l ‘ 

-.00179 

-.00160 

-.00143 

-.0012 9 

-.00117 

-.00106 

8 

-.00273 

-.00237 

-.00209 

-.00186 

-.00165 

-.00148 

-.00133 

-.00121 ’ 

-.00110 

^.00101 . 

10 

-.00242 

-.00214 

-.00190 

-.00170 

-.00152 

-.00137 

-.00124 

-.00113 

-.00103 

-.0CG95 

15 

-.00180 

-.00164 

-.00148 

-.00164 

-.00123 

-.00113 

-.00103 

-.00095 

! -.00088 

-.00081 

20 

-.00135 

-.00125 

-.00116 

-.00106 

-.00099 

-.00092 

-.00086 ^ 

< 

-.00079 

-.00074 

! -.00069 
1 . 1 
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TABLE AIV.- AUXILIARY a^NCTIOK Fi + F2 FOR DETERMINING 
DCViflAmGH FJilCTION F AT IffiGATIVE V/iLUSS OF x* 

. . [F(^x^ = F, + Fg F(+x*) ] 


y* 

F 4 . F 

y* ] 

a, F 

‘ .T ^ 2 , 

0 . 

4 . 0000 

24 i 

t 

~ 0 , C070 

2 

-1.3333 

! 

26 I 

-.0059 

4 

-.2667 

i 

28 i 

-.0051 

6 

-.1143 

■ ■ 30 , 

-.0044 

8 

-.0635 

32 

-.0059 

10 ! 

-.0404 

1 

. 34 . 
1 ‘ 

' -.0035 

12 

j -.02 80 

13 

-.02381 

14 

•-. 02.05 

37 

-.00292 

16 

-.0157 

1 

61 

-.00108 

18 

-.0124 

93 

1 

-.00046 j 

20 

-.0100 

117 

1 ’-.0002 9 j 

22 

-.0083 

141 

j -.^00020 I 
! i 
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APPEWDIX B.- PEHTIWEWT IWECRf-IA.TIOW AWD COMPUTIWG 
PSOCSDURES res Uffii WxTH THE MUTTEP.FEKL MTHOD 


A lifting line used to represent a w3ng ia pLaced in a p^osi— 
tion coiTo spending to the quarter-chord lino of the wing, '[he 
distribution of circulation along the lifting line is exp:rossed hj 
equation (2) . 

Wo generalization can be made as to’ the n-jnbor of teras which 
must bo retained in the series to ensure acceptable accuracy. 
Mutterperl implies that four are sufficient and utilises this number 
in <all applications. It should be noted that all loadings predicted 
by the series as it stands will be symiO tried. In addition, 
equation (2) cannot satisfactorily approximate a curve containing 
discontinuities such as would bo prod viced by flaps or ailerons. 
Mutterperl mc,de no comment .as to additions or alterations to the 
series which would enable circumvention of these limitations. As 
.a result, while it is believed that such modifications could bo 
included, it is not known to wh?,t extent they would Increase the 
complexity of the mathematical evaluation. 

Since it con bo shown that in a theoretical appx’O.ach using a 
lifting line at the quarter— chord line, the d.ownwash angle at the 
three— qiiarter— chord lino most closely a.pproximate s the true angle 
of attack of the wing, the control points wore placed a.long this 
lino. The number of points required is dictated by the number of 
unlaiown coefficients retained in equation (3). The location of 
these points spanwiae is app;arently ai'bitraiy; however, since 
Mutterperl placed them on thu right wing half at r, = 0.1 74, O.pOO, 
0. 766, 0.940 ( = 80°j 60°, 40°, 20°), this arrangement has been 
followed in all applications presented herein. 


Deter:aj.nr.tion of Downwash 

The expression for dovmwash at a control, point, as determined 
thr'ough tho Blot— Savart relauion, is 


3 ^ 
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sin a 

n=0 


Ar 


ii- 


_ ) (2n+l)a2n+i |2n+l V 3 r BL./^L+S'^Ly 

Tt sin (gn-g) '/ + _sl r^^ cos fgn+l^-^ cos cpd^ 


s’ cos A sin t 


®1?Jo 


”1 ^ A"^ 

— i -^ + tan A 


Be \B3 




,•0(0 


cos (2n+lyCpd:p _1 ' ^~B' . o i- tan A ; 

\ 


V«-’r26 


" ! - V® .it 


nn/?co3 (2n+l)cp‘ cos cp-— ^ 

• c^r \ 1 ' ' \ Si 

F(C+F) 


dCD 


I 


. £i_ ir*i - ter. a'i i“'‘ 

B2-./- G Bl\BL A/o ^ 


+ i- 

s’ V 


(AT-y') tan A\ cos (2n+l)cp<3-? l 

1- 11 ^ rJ_ P_ 


(Bl) 


cos cp + 


Equation (Bl) may Be reducx^d to a simple expression, cmuau.n 
ing only the unlmoiras a^n+i ^nd tlieir n'ornerical coeff i^^-ents^ 
simply By the introduction of wing geometry and the geometry of^a 
control point. Such a reduction sho-old therefore Be carried cut, 
at each of the points. 


Solution for Additional Loading 

Since for such a solution the wing is considered a flat plats, 
all geometric angles of attack Become sin a and the ccnsuanu 
factors in equation (El) Become one. Simultaneous solution of 
these equations then evaluates the unknowns agn+i^ whicn are ^ 
suBseauently introduced into equation (2) to produce tno aesirea 
expression for additional loading. The unknowns can also Be intro- 
duced into the following ex?pressions derived from equatxon {dj, to 
obtain the values indicated. 
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cICt 2jr2bC-, 4ir2ai 

Q.t i= 

da 3 1-rX 


(B2) 


£i 

Cl 


?(!+>>- ) 


rta.i[ 1— (l— a) cos cp ] 




sin 9 + as 


cicr. 

Cl cav 


sin cp+ag sin 39+3s si'o 


T'icp = 


ii_ ( 1 + 3§a. _ 25_ 

3jt \ 5ai 7ai 


+ _az_ 

l^a-i 


sin 3 ?'^a 5 sin 


^cp+Oy sin 7?) 

) 


P^+ct' sin 7(p) 
(S3) 
(b4) 

(35) 


SclutiOii for 3a sic Loading 

As in Fallnier, the hasic leading on a twisted and/or cambered 
airfoil, can best be calculated by detennining the total end addi- 
tional loading at some finite lift coef.ficient end subtracting the 
latter from the former. Tlie procedures involved are parallel to 
those of the Falkner method as well. An arbitrary angle cf a.ttack 
ag can be selected for the root section of the wing., from which 
all local angles of attack can he measured. Svlstitution of 
sin ctiocal f^r w/V and sin ag for sin a in the expression (B1) 
will ros’olt in a set of equations which may he solved simultaneously 
for the val’ies of the coefficients a^-n+i . If and sin Os 

are then introduced inoo the following, the coiroct value of the 
lift coefficient for the wing at the descrihed attitude a will 
result 


Cl 


_ 1 sin a 

1+X 


(B6) 


In addition, if the values of sin ag are substituted 

into tile following, an equation for the ordinates cf the cui’ve of 
total loading on the wing at a results. 


cicri 


4SCl 

jtaih 


(hi 


sin 9+3.3 sin 3<P+Ss sin 59+a-7 sih 7®) 


(B7) 
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If the values of coefficients ottcined for a solution for the addi- 


ordinates of this last from those of the total loading curve will 
give those of the desiirod oasic loading cxurve. 


As in the mothod of Falkner, the error introduced into the 
solution hy the ass+miption that all section liftr-c\jrve slopes 
were 2« can be readily eliminated. If the actual cj^ does 
not vary across the wing^ the coefficients a^n-i-i should bo muxti— 
plied by the ratio of actual to 2jt. If the actual value doss 

vary along the span^ this ratio should be included in equation (3) 
as a function of the spanwise variable 9 . 


These instructions apply only to unyawed wings devoid of camber 
and/or twist. 

If the local angle of attack sin a is introduced into 
expression (Bl) in place of the dorawash ratio w/V, this expres- 
sion can bo wrcitton 



Correction for ci^ 


Computing Instructions 



n=0 

where 



2n+l Bp vui2p-!-B'T\ 


bl 


7l s in ( 2n-i-l ) ^If 

s' COB A sin -ijr 


(B9) 
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Computing fora B(1) is used to calculate all factors which remain 
constant thrcug^iOut the summations for any one control point: Ag, 

Al, Bl, s', C, Fg, F4, Fg, Fg, Fy , Fg. Tlie computation fora 
for those factors, which vary with n only throughout the summations 
for any one control point, is presented in form B(2) . Computing form 
B(3) is used to apply Simpson's rule to the integration of Fa', F/. ', 
^5 S ^"7' ^nd Fg ' . Hie factor Ksn+ij which is independent of. 

wing shape and so can he a.pplied to all wings, was calculated from the 
relation 

Wl (310) 

Jq cos Cp — cos ly 

and is pi-esented in the following table : 


j Control ' 
! point 

-j.- 

, J ■■ ■ , 

1 

n 

i 

r~r-| 

Ksn+i i 
> 

I 1 

i 1 

0 

2.54928 i 

i 

1 

1 

8.26^92 ! 

1 

j 

2 

8.87752 j 

! 

3 

6.10976 I 

i 2 

0 

2.47976 1 

i 

1 j 

1 4 . 5S64 o j 

1 


-1.90944 1 

1 

; 3 1 

i -4.63656 ! 

■3 1 

0 

1 2.33120 ! 


I 1 

0.47920 1 


i 2 

-3.47696 ■ 

j 

i 3 

3.39776 1 

h i 

i 0 

2.00112 i 


1 1 

-1,92340 j 

1 

! 

! 2 i 

1 1.38068 i 

* i 

i 3 ! 

i i 

! -0.53112 

1 


The results of the integration are presented in form 3(3). 
Form B(4) is the form in which the components F are multiplied 
together and the results are summed producing fom* equations, 
one for each control point. La.stly, the form for the sim’oltaneous 
solution of these equations (reference 9) and the resulting values 
of the unknowns are presented in fora B(5) • 
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APPENDIX C.- DEVELOPMENT AND COMPUTIK PROCEDURES 
FOR USE WITH THE ^/ffilSoIWGER METHOD 


In the German reports a.yailahle, which described this method 
the mathematical development was not complete . It was thought 
advisable therefore to include the development in the pi-esent paper. 

From lifting— line theoi"y the equation of downwash at the one- 
quarter . chord line of a straight wing is given by 


This integral eqxxation is solved by Multhopp by an integration 
foimula (reference 5 ) • 

The equation of downwash at any point xy of a straight wing 
is given by 


Weissenger divides this integral into two integrals j one of which 
is the same as equation (Cl) which he solves by Multhopp 's method; 
and the other which he solves by a method analagous to that of 
Multhopp . 

Tile mathematical development is as follows : 

With G = r /bV, T) = 2y/b, q = 2y/bj ar = b/c;^ 

and setting x equal to the distance to the three— quarter chord 
line X = Cq/2, the equation (C2) becomes 



(Cl) 




4it 
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or 


V 

V ■ 


ri 1 



g-LCiOln 

7-7 



La [ar(ri-ri)] G'ti d'P, 


(C3) 


where 


La= 0 [cr(n-r))] 


*/l+f ar)^(ri— — 1 
ar(7-ri') 


(C 4 ) 


or 


La= 0 (v,!.i) = L[ar(cos cp-., - cos cp,.)] = 

' arv(cos cp^ — cos cpp.) 


of CO, 


TliG first Intogral of equation (C3) can be vritten as a function 


2it ^Li Q— q 2tc Jq cos cp — cos cp-^ 


(C5) 


where 


cos cpy = cos 


vn 

m+i 




An integration formula, gives 


^ f (7)dq = V f (lin) sin cp^ 


y_x 


m+l / 
n^ 


(C 6 ) 


whore cp^ = L(qn) is the value of f(ri) at r\^r^. 
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Equation (C6) holds exactly if f(q) can he represented hy 


2m-l 


2ra 


V=0 V=1 


f(q) = 2^ Avt)'' = > av sin vcp 


or, in addition. 


fC'-p) = 

n=l t^x=l 


Letting G(cp) = f(q>) 
then 


m m 


G(cp) = ^ Gn 2 


n=-l |-Li=l 


(CT) 


^ ) G„ ) di sin dicrvix cos dicp 

dcp m+1 i_, ••' l_ 

n-1 di=l 


and the integral in equation (C5) beccmes 


m 


m 


) Gn r ' 1^1 C03 ji gaL. dcp 

m+1 L • Jo Z_. ‘ 


/ 

n=l 


di=i 


cos cp — cos 4>v 


Now 


r” dcp = ^ -'^^0 

Jo cos cp -g. cos cPq sin cp^ 


(C8) 


or 


r V (-^1 Bin mcp.) f £ P? - Hicp ,. . ^ dcp = n V s . in ,. ^ j. cpT) sin . ^X i Cg . ^ . 

J(o Ls ~ Vcos 9 - cos 9v/ /_j sin cpv 


(1 1=1 


L ) 

Pl=l 
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30 that equation (C|5) hecomea 


(!') t 

ri 

m 

\ 


sin M-iCpn sin l-iiC[)y 

'■siv' V+i/ L 




sin 9y 

m=l 


di=i 



Fi 


in 



- ^ 


Y 


sin PiCp^ sin PiCpy 

m+l L 

L 


sin cpy 


n=l M-i==l 


(C9) 


The |ii aeries is independent of wing geometry and may he put into 
a. coefficient usable for all wings . 


From equation (C9), Tor n=v, let 


bv,v = 


m 

—A V 

(lii+l) sin <Pv /_j 
Pi=l 


uisin^ PiPy 


and for n ^ v, let 






-1 

(m+l) sin cpv 


I 


Pi=l 


Pisin PiCpn sin PiCpy 


Then equation (C5) becomes 


(CIO) 


(Oil) 




(Note: The summation prime indicates that the term of n = v 
should not be included.) 

Expression (C12) gives the induced angle of attack on the one— 
quater-chord line, at the span station v, in terms of the siaamation 
of n spanwise values of the dimensionless circulation. 


) ^ ,n 

L. 

n=l 


(C12) 
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Equations (CIO) and (Cll) may "be simplified- with the aid of 
the summation formula 


m , m 

) l-ii cos M-iX = real, part of( ) • 
Z-j ' \Z_i 


_ m cos('m-fl)x - (m+l)c os mx + 1 
2(cos X— l) 


or 


bv^v = 


ra+l 

4 sin CPv 


(C13) 


t)v . n 


9 n I 

(cos cpp - cos 


l-(-l)^^ 
2(m-:-l) - 


(C14) 


The second integi-al of the integral equation (C3) ex' 



■LA(Vj!-i)G' (rj)drj 


(C15) 


is solved in an analogous fashion. 
The integi’ation fonnula. gives 


I f(q))dcp 
Jo 


= .JL. r f(To)+frTi'4-:iT 


M+1 




u=l 


where 


% 


■QJL. 

M+1 


(C16) 


f 
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Dimensionless circulation is given by 


m m 



n=l Ux=l 


sin sin Mi© 


or 


Letting 



2 

m +1 



n=l 


/ 

Z_. 

1 ^ 1=1 


M-xsin 




9jj 0 03 1-1x9 


f J9) = ' ) Hx-sin P-iiPn cos 11 x9 (CI7) 

n+i /_, 

lii=.l 


then 


^ Gn ?n( 9 ) 

L~i 

n~l 


Expression (CI5) may liiodified so that 

-L Ly\(v,a) G’(7)ti'n = - I L^(v_,u)G'(©)d9 

2 « J_ X ■ 


m 


- ) 

L I 

n=l 


i'a J 


nTt 

: 


Then applying equa.tion (Cl6), ^ J La(''’,ii) G’(rj)dr| 


L/^(v,|.')Gri fn(9)'| d9 

J 

n 


n 


M 


1 

J- 

2 n 


n=l 


G- j fl'AC V , 0 ) .K f n, o.)]La( V,_ M+1 _ ) 
M+1 ^ C 2 


+ ) r 

-L_ ^ 

n=l 

- Y Gn — -j— I , q) + ^La( v.M+llllfj^-taJ_. ^ ^ LA,(v/Ofn p 1 
L 2 (M+ 1 ) L .2 ‘ J 

n=l 
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For simplicity let fn,!--. = where 

equation (C17) for cp = = pn/M+1 


rom 


Then if 


„ f [LA(v,o) ](fn.o)+[L A(v>M+l)](fn.M+i) + LA(v,n)f„ u 

SV^ - 2(M+1) \ 2 L. ’ 

p=l 

(CIS) 


expression (CI 5 ) becomes 


m 


y Gngv^n 


(CI 9 ) 


n=-L 


The solution of equation (C3) is times equation (C12) plus 
arv times equation (CI 9 ) or 


w ^ 2 ( 
V A 


- 2 Gy + 



1 


\ 

m 

\ 

■ L 


Go. 

) 

y 

+ ) arv Bv yTL 

II 




n=i 






ary 

Sv,v 

A 

+ 

y a-v Sy^n Gn - 2 \ b 





li 

II 


v^n'^n 


so that 


m 


^ = (2bv.v + arv gv^v) Gv — ^ (Sby^n— ary Sv^n^ *^11 
T /L. 

n=l 


^ = b'''v Gv> 


f ' 

; b*v,n Gn 


/. 

n=l 


(C20) 
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vhoi-e 


b*v = 2bv,v -i- £5a-v gv,v 
^'"v,r! - 2b-\^n~ Sv,n 


(C21) 


Tlio prime on the summation, indicates the.t the .term of n = v 
should not he a.6.ded into it: ' 

Equa.tins the doT-mwasli to tirie local a.nglc of attack of the plate. 


Cv - h-y “2^ ^'^‘v,n On 


(C22) 


For a swept wing the equation of do^mwash at eny chordwise 
point for y>0 is given by 


1 pb/fc ri(,) r 

* /II rt?) 






7 - If I A)^+ iy-y)^ J 

2L^-iL^n,A_ dy 


<iy 


i-b/s [(x — h7i ta.n A}^+ (y— y)2]3/2 

h/s 

’^^Jo 


X - y ta.n A 

[(X - iyPS. A)S, (y-y)5]V5 'll 


(C23) 


The f'irs'b integral is the dovnwa.sh due to the trailing vortex 
sheet and the last two integrals represent the velocity induced 
hy the liTting line. By integrating the last two integrals by 
parts a.nd reorra.nging into dimensionless quantitioSj the preceding 
integrals may be put in a form similar to equation (C3) . 


a = 


2:t 


nl ^ 


r— 'q ' 


or 


pi 


2^ J^2 


'^/Sn,T)) 0'(q) d 


(C24) 


(Woto that in. tho^f ollowing oqua.tion (C25) the squaros under the second 
radical of LyXilyh) fo-" aro summod. In roforonco 3 these 

wore erroneously sho\-n a.s a difference . ) 


where foi- t)<0 


LA(h^^) = — ; t/fl..! ..ar. tp,n. A (V | t) | ) ] ai- ) ,T]-ri ) j _ J., ^ g. tar Ar -;- 1 

ar(ri--Q) t 1 + 2 t) (ar) tan A J - - , > 


1 + 2',) (ar) tan A 


and for r[ >0 


ar (ii-vf) 


)(C25) 


Equation (C24) with LA(l^^i'|) ^iven hj'' equation (025) io for y>0 and will give 
vclueo for only the positive span stations. For values with y<0, expressions similar to 
equations (C23) and (C2^) aaust he derived. Values of L^(T)jiy) ere needed for y<0. if an 
unsymmetrical wing is to analysed, “■ 

Equation (C25) with q = cos cpy and iT = cos Cfp, and ai” at span station v as ary 
fox’ cos cpy >0 hecomes, for cos cp|j, <0, 




ary (cos cpy — cos CpT 


Jj fen A (coa, cpy — jcos cptj. !)]"’+ (ary)^(cos cpy - cos fqj)‘ 

1_ 1-2 ary cos cpy tan A 


-xl 

J 


+ ^ [1+ary tan A co o cp^]^ + (ary)" " cos^ 

1+2 (ary) cos cpy tan A 


'>(C26) 


and for cos cp^ >0 


La(v,p) = V[l+aJ^v tan (cos cpy - [cos cpj, [)]^ + (ary)^'(cos cpy - cos cf^)^ -1 

ary .{cog cpy — cos 




H 
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Equations fox* determining Gy at spanvise points of 
^ =■ cos for any ving may therefore be detemined from 

equation (C20) . 

Summarizing the computations, the relation equa-ting downwash 
to the local slope of the plate at m points along the span giving 
m equations with m unknowns, G^ is 



where 


ov = angle of attack at span station v 



b*v = 2bv^v + Sv,v 



b 


m+1 


4 Bin cpy 



Sv,v = ev,n « = V 


2(M+1) . 


i 


La( > o) ( fn . o ) '• ^a( ,M+1 ) 


2 



Jl = winr; span 


•cy chord at span station v 


La(v,m.), for a straight wing, see equation (ci*-) . 
La(v,|i), for a swept wing, see eqmtior. (C2c) . 


f 
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m 

^n,n = ^ si’’ I’l ^n cos Hi 9^ 

/. 

M.l=l 

where agcin 


9n 


-JDiL 

m-fl 


Cf>v/ 


VJt 

m+1 


Cp^i 


tiTt 

M+1 



The coraputations required in the preceding group of equations 
may he simplified if such values as hy v n various 
and fn^p various m^s and M*s are tabulated. Then o so3.ution 

for any wing consists of a substitution of wing geoiietry Into the 
L^(v,n) function, equation (C^) or equation {Q 26 ) , and a coatination 
of the tabulated coefficients tc obtain ra simultaneous equations 
with m unlcno->TO3 G^. 


The computations for a. s^/JicDistrically loaded wing may be still 
furtner reduced by an alteration to the preceding equations and 
coefficients. For a sy.nmetrical \7ing with cr without camber and 
twist, the distribution of local angle of attach is synmietrical 
about the center of the wing or 


then 


and 


^■v + Oto-i-l-V 




o-v = 


V 


goes from 1 to 


m-fl 

,2 


n goes from 1 to m 


H goes from 1 to M 


‘5ij- 
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The s’jmination terms ca.n he vritten siich that n from 1 to m 
hecomes n + m^-l— n from 1 to and |i from 1 to M becomes 

M. + M+l-ii from 1 to (For the cases n = and \x = 

the coefficients are expressed as their former values as will he 
seen directly.) 

Now if equation (C20) is expressed as 

a+i 

'(L 


Or. 


V - ) (^'V,n ^"’V,m+i— rJ 

n=l 


n 


\h',eve 


n ^ S±1 


V • 1,2, 


m+1 


Then 


where 


= B*v Gv 


ra+1 

2 

y ^n 

n=l 


(C28) 


B\',n = + h-\,,+i_^ V = 1,2 . . . n 4 


(C29) 


(C3C) 


These coefficients may also he erpressed a.s 

3'^; = 2hv,v + £’-'V BV,V 
B'*^v,n = ■“ Sv,n 


where is limited in that + ^,m+i-n ®v,n = ^,n ^ = ^ 


To find g consider the expi*ession for g 

v,n ^ ^,n 


8. 


2(M.l) 


[l-A(v,o)](fj^ q) + [l‘A(VjM+l)](fn j4+i) SP T /• w ■,^ ^ 
+ ^ L^(v,M+l) f^^ 

" U=1 


The summation term of g^ ^ can be written for n f 


, m+1 


14-1 

2 


M 2 

Z' ! " V ! n.u'''^A n,M+i-nl 


H=1 


d=l 


[‘ 


" 2 




for n = 


m+1 


Me l_ ^ 

Z! ^rx,p = + Ia(v,M+ 1 -p) f^^M+xJ + (C3I) 

d-1 U=1 L 2 J 


M-1 

2 


VJ1 

U1 
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vhere as "before 


m 


m+T fe+T 

|jii=l 

examination of equation (C32) will indicate that for n = 




= 0 


for n 4 

2 


^a( i<A\ ^>^^2 j ^ai+i— n,^^^ 

\ /' 2 \ / , ^ 

have values onl.y- for .even Ui, tut for even P-i, fn,n = " ^m+i— n,iJ. 
so that these terms equal zero. 

Then equation (C3l) becomes, for n 4 

M ~ 

^ iA(v,|.i) fn,p ~ ^ ■'■ ^m+i— n,M^ La(v,u) 

^ 1=1 (. 1=1 

+ (^n,MH-i-y. + fm-M— n,M+i-(J.) La(v,M+1-u) ] 


and for n = 
M 


x^l 
2 
V 


V- 


1-^ 


La(v,h) fn^iJL = ) C(fn,li) La(v,u) + (fn,M+i-ia) lA(v,M+l-ix)] (C33) 


Further examination of equation (C32) will show that for only the 
odd (ii 


i’n,(J. = ~ fn,m+i-u 


However, the second summation in equation (C33) has n = in which instance the even 

terms in equation (C32) vanish, thus the second summation becomes 


M-1 

2 



h=l 




La(v,M+1-u)J 


The first summation in equation (C33) has terms of fn u = " ^n+l-n u M+l-u = " fm+l-n,M+l-u 

for even or the coefficients of La(v,h) ana Ly^(v,M+l-^l) vanish for even 

In addition fn,u+^m 4 -i-n,u = ‘ ^n,M+l-u-^^m+l-n,M+l-u summation is 

M-1 

2 

\ (fn,u+m+l-n,u) [hA(v,h) - La(v,M+1-h)J 

where fn,ia‘'‘^M+l-n,|i is obtained from (C32) for the odd terms of 
Lastly, 

(fn,o)LA(v,o) + (fn.M+i)LA(v>M-*-L) _ (f'n.o+i‘m+l-n,o)hA( v,o) + (fn,M-‘i +fm+i -n,M+i )La('^,M+ 1) ^ ^ m+1 

2 2 ^ ^ 2 


= ^^n.o~*-Wi-n.o) [la(v,o)-La(v,M->-1)] 

2 


for n = 


m+1 

T" 


where fn,o-^f'm+i-n,o is taken for only odd Ui . And, for even o = i'm+i-n o' 

a factor may be expressed as * ^ 

^n,u = ^n,u for n = ^ and odd 


So that 


m+1 


fn,u = fn,u+fm+i-n,u for n / and odd 


VJl 
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li=0 a.nd odd |Jii 


f or n 

f" ,, = £hji.£. for n = p.=0 and odd Hi 

“j. o o 


and 


M-1 
2 

[La(v,u) - LA(v,M+l-tO] 


’ 2(M+1) /, 

lJ=b 


Examining equation (C32) 


f ,, ~sin lii 

iti+1 


fm+i-n,!i ~ sin m 


sin fuii« 


u+1/ 


liJjiZL _ sin 

m+l m-'rl 


= sin M. 11 T cos — sin sss Ui« 


or when (ii is odd. 


Thus 


„ . ij. 1 n« 

Am+i-n,(i ~ sin -— 


■^n,!i - -^Tn+i -n,n 


Tlien, as before 


M-1 

P 


GV,n = 


-1 


2(M+1) 


f n,|i [ La( V , d) - La( V,M+l-d ) ] 


d=0 




(C3^) 


9i+t'C * ild' VOITSI 
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where finally 


iriji.!. = n = S|1 and odd Hi 

fn^lJ. = 2fn^!i for n ^ and odd hCi 

fn,ia = fn^ii for n f !^ = 0, and odd p.i 


(C35) 


fn,|i = for n = ,u = 0_, and odd |ix 


and 


m 

’n,!i = y III Sin di ^ cos y~ 
■* m+1 i_, m+i M+1 

Pi=i 


A further simplification to Sv,n also he effected. From 

the hlnomial theorem; 


^ 2 a. 8 16 aA> 128 


. a r 1 + 1 'i'f - 1 (kV + 1 

! 2 .a/ o \a/ lo 






12H V,v 


for a>b 


(C36) 


and 


. y ,,2 . - 1-1 I = S '- J. a ' X - - 

-vd. -hD _-b.! + 


1 a.2 1 a'i ,1a® 


= a 


i + 1 


1_^ 

2 Vo /a/ 


_ 1 /0_'\ 


1 ' 1 




8 Vb/a 9 16 '\bJ&J 128 Vb/a 7 


’or a <b 


(C3T) 


Wow letting 


T(t) = V 


- A.^ a^ — b ~ 


6o 
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then for t = ^ 

O- 


T(t) = .1 + ^ 




" iIj 


t® for t <1 


or 


n-ft') - t + - i ' -1- ^ 

nt) - t + ^ ^28 ^.t.- 

Equation (C3^-!') ^^a.y he ■V'rrittsn 




M-1 

2 

8v,q " 2(M+i) / , 

n=b 


(039) 


wnere 


^a(’'-^^ 1^) = ^A('h'-i) “ Ia(v,M+1--+0 
This can also he T-7ritten 

ALa(v,u) = LA(q,h) - LA(3.-h) 
so that from equation (C25) 


:C4C) 


(C4l) 


ALa(v,|j) = 


A-'^r I-t-ary ( n-'n ) tan A ] [ ary ( q- q ) 3 




arv(-'i-A) 


1 /~i 1+ary ( q-l) ta n A j ^ + [ ar., . ( n+^ 

l+23ryq tan A aA_, (q-i-’',) 


tan A 


l+2arvq tan A 7 [l+arvq tan A]' + [arvq 


-.2 
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Then 


ALa(v,h) = T(ti)- 


— -I T(t2) + T(t3) tan A 

l+ 2 arv’nta.n A j 


- T( tg ) tan* Ar- -i- : -~ 25 - 

ai’v vt) 2 — / 


(c43) 


vnere 


tx = 

tg = 

ta = 


arv(ii-T]) 


-r tan A 


ai\)(T)+nT t+T] 


tan A 


ary r; 


i+ary T) tan A 


and T(tx) takes the sign of (ii— ?i) 

T(t3/ takes the sig-i of ta 

Several of the functions cf equation (C 43 ) are independent of wing 
parametei’s and ma.y be tabialated for varions q and 

8 v = Ke=l, 


K,. = 


ri=— r]« 


Id-, K = am 

Traf ^5 _ .■=■ 


ri+T» 


where, as befoi'e. 
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Then 


AL/^(v,m.) = T(ti) 


where 


1+2 ten A 


[ T(t2) + T(t3) tan A] 




- TCtg) tan A - Py ^4 


”^1 = |Pv + tan A I 


4 - _ 

^2 = 


j p y Kn + ICg tan A | 


and 


Py Kg + tan A 


T(ti) tahes the sign cf Kj 

T(t 3 ; tai:es the sign of -Cg 

In sammaryj for the syramotrically loaded wingj 
tn+1 
2_ 

n, V = 1,2 . .. 3±1. 


ay 

= E*y Gy — ^ 


H=1 

B-y 

= 2hy^v + arv gy^y 

B%,n 

“ ^^v,n “ Sy^n 


f By,n By^ja+i--n 

Bv,n 

Ibv.n fcr n = 


7 mi 


■'V^V 


'v.n 


jTi-:-l 


4 sin cpy 

sin cfa 


(cos cp-g — cos cpy)^' J I 


i [ Mzli 
2m+l 


a— y -| 


(C44) 
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Sv n taken from equation (C39); is taken from (C44) 

an^ fn,(i is from equation (C35)« 

Limitation of the Series 

The nurn’oer of coefficients Gn required for accuracy depends 
upon how rapidly the series equation (3) converges. Weissinger 
used m equal to 7, I 5 , and 3 I in his investigation and concluded 
that the results ohtained with m eqToal to 7 were nearly as 
accurate a.s those with m equal to I 5 or 3d. For this reason m 
equal to 7j or four coefficients have heen used in all of the applica- 
tions of this method presented hei’ein. Tlie niamber of tenns required 
in the interpolation function fn,u must also be established. Again 
Weissinger used M equal to 7> 15 ^ s-^d 31 snd foijnd tjf^t results 
with M equal to 7 proved as satisfactory as those with M equal 
to 15 or 31 * Lastly, it should be noted that equation (3) cannot 
satisfactorily approximate a curve containing discontinuities; 
however, a modification which will enable it to do so has been 
developed by Multhopp (reference 5) • 


Solution for Additional Loading 

Since in a solution for additional loading the wing is 
considei-ed a flat plate and all angles of attack av are equal 
to a, equation (C28) may be modified to 

ra +1 

'2 

(C45) 

n=l 

Evaluation of this equation at the several stations V produces a 
set of equations containing the unknown circulations Gn which can 
then he solved simultaneously to obtain the values of these 
circulations. 

Substitution of tho values so obtained into the following 
expressions^ results in tho values indicated: 

m 

\ Gn sin cpn (C46) 

da m+1 Lj a 
n=l 
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or for a symmetric wing, 
Gi 


m— 1 
“2 ‘ 


da m+1 \ a, Z_) 


^cp = 





n=l 

- 

O Q *V% < 


1 

CL 


a 

dCL,/da 

ClCv 

- OAT? 

Gv 

1 

CLcav 

a 

dCL/da 

iference 6 for 

It 

B 

0. 

3:32i+Gi 

+ 0 

. 903 G 2 • 

0. 

3827 G 1 

+ 0 

. 7071 G 2 


(Cl+T) 


(Cl+8) 


(ci+9) 


(C5C) 


Solution for Basic Loading 

The hasic loading on a wing with camber and/or twist can be 
determined in a manner exactly paral3.ei to those of Wallcner and 
Mutterperl. An aibitr 8 Jr 5 ' angle a^ is sexected for the root 
section and the values aiocal^ raeasxired from it. If these values 
are then substituted in eq.uation ( 028 ) , and ix the resultuig 
equations are solved simultaneously, values of Gp. will be 
obtained, which when inserted into the iolxowing expi'ession wxl^. 
give the correct lift coefficient for the wing at this attitude. 

rn— 1 

ot = ^1+1 + 2 Y On sin -Dr, ’ (C-'l) 

m+l\ 2 z_.. 'V 

n=l 

If these values of Gp are also substituced into the following 
equation, an expression for the total loading cur\'e will result. 

(O'; 2) 


01 - 2 \) Gv 
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Of tho valuoQ of Gn obteincd frora cc,lculations for tho addi- 
tional loading aro then eubetitutod into equation (C52), ordinates 
of the additional loading curve will ho obtained, which when 
subtracted fron thoso of the total loading curve will give those 
of the dosii'od basic loa,ding ctirvo. 


Correction of ci^ 


As in tho other methods, the eri-or introduced bj’’ the assumption 
that tho section lif't-curvo slope is in all instarices 2jt can be 
roadily corrected. Ihe correction is accoraplishod by modifying the 

specific circiilation ordinates Gn by the ratio ^ ^ wher’o 

2n 

specific values of this function a. t span s to. t ion n must bo 
dotemined if this function varies along tlio span. 


Computing Instemctions 

Tho functions are dcte?miinod on form C(l) for a 

swept wing and fozm C(2) for a stral^t wing (Ac. 2,0c = ^) • 
both cases tho coefficients Ki, IC2, Ka, IC-., and K5 aro obtained 
from tabic Cl. In tho values of T(ti), T(tn), TCt'a) are obtained 
by entering chai't Cl with the values cf tj., ts frem columns" 
(10), (11), and (?-2) foim-CCl). Form C( 3 ) contains those computa- 
tions which result in tJ.io p3cprGS3ion3 containing tho 

unknowns Gj^. Tho values gv.,n in •column 9 of form C( 3 ) oro obtained 
as follows: Consider the values in col^mm 4 of this form a.s four 

groups of f otir membors with the grovips identified as 

A ■’.-.’iion n = 1 

B i-dion n = 2 

C when n = 3 

D when n = 4 

Similarly, tho values in col'jmn 5 of this form can be considered 
grouped a.s 
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1 when V = 1 

2 when v = 2 

3 when V = 3 

4 when v = 4 

then 

= r A X 1 

51 . 2 = 4. B X 1 

gi^3 = E C X 1 
Si,4 = K X 1 
Ss^i - ^ X 2 

5 2.2 = ?' B X 2 , etc. 

Ihe conpirtrijig form fon the siinultaneous solution (neierence 4) 
ef the equation (C4-;) is given in form C (4) . The equations are set 
up as follows: The first equation consists of the first four nuabers 

in column Ip foim C(3), the second equation the second four, and so 
forth; the first number in each group being the coefficient of 
Crx/rt, the second being the coefficient of Ga/a and so forth. 
Simultaneous solution gives the g/cc's with the corresponding spflfl. 
station. 
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-A- = - 45° — -/.007 

X = ‘^7^ I- X 

M - K = = Z.OS7I. 



z 

J 

4 

S'' 

6 

1 

a* 

3 

10 

II 

12 


14 

15 

7? 

17 

18 

19"" 

20 



/- 

(?-A) 

V 

A" 

wa; 


(4)Kj 


ianA 


1(6) 1 
1 -f<9)\ 

4= 

T(i) 

7(4) 


tanA. 

/ 

/f2^ 

W 

m-i- 

(/?h 

= 0^ 
-(3^ 

-0(n 

-or) 

/ 

o 


.205? 


jono 

Ji2Z9 

-2.6/23 

.039? 

57/2^ 
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/.a/ 
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Form C-- Computing form for Weissinger's 


METHOD 


(Underscored numbers are sample catculcrtions.) 



"2 

3 

4 

5 

6 

7 I S 

9 

lo' 

1 1 

12 

13 

14 

15 


n 

/’ 

-fn/A 

Taslb- 

CJL 

iUi^) 

(zol, 

or 

(io\^ 



n 


/i>^~ 

(9>6o) 

2 6pn 

TaqU 

cm 

(12) 

-00 

ViiSLC 

cm 


1 

1 

O 


'2.330 


I 

1 


-,Sa?6 


0 

-J0^4 

!0.4SVi 




1 

-1,414- 




z 




33Z84 


0 




z 





3 




o 


_J2_ 




3 

,4l4 



4 




.i^za 


C? 


7 

2 

0 

-A4j4 



2 

1 




Z.O7Z0 


Z7 




^ 1 





z 




0 








3 




^ 



4 





n 

Gi2c^Qi 

(K (X 

^PflN 

jrATiaS 

COS-^ 

/ 


.4254 

2 


.70 71 

3 

. 54^4 

.7327 

4 

,U7^ 

0 


1 Equation 

1 

2 

3 

4 

CoNSTAMT 

L 


L 

1 

'^i/oc 

WFm 




j Gz /ot 

rZm7\ 




1 ^j/j< 

.?7Z0 





-,!3CZ 
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TABFE Cl.- CCF8TAIIT FACTOE PSCilTiEED FOR TEE 
CALCULATION OF ALa(Vj(Jl) FOE n = M = 7 


V 

1 

Ki i 

Kp 

K3 _L 

K4 j 

i 

1 /n-- i 

i/n+iT 

.. i/n 1 

2Tf/n^~2 1 

Ti-i/n+i'i 

*1 

— 

o| 

-l^.i 4 o 6 1 

0.5198 

1.0824 j 

-13.6610 1 

-0.0396 


1 

I 

c» ! 

j 

.5412 

1.0824 j 

00 j 

1 

0 


2 

4.6125 i 

.6131 1 

: 1.0824 I 

3.9992 I 

.1329 


3 

1.8477 i 

1 

.7553 ' 

1.0824 ■ 

1,0823 1 

.4142 

2 

— t 
0 , 

- 3 . 4 i 4 i 

.5858 

1.4142 1 

- 4 . 0000 ! 

i 

-.1716 i 

1 

1 

1 

1 

-4.6125 

.6131 

1 . 4 i 42 j 

-5.2254 i 

-.1329 I 

1 

2 

00 j 

, .1071 

1 . 4 i 42 i 

{ 

0 : 


?! 

1 

3.0S26 

.9176 

i 1.4142 1 

j -1 

i 2.1650 ; 

\ i 

.2977 

3 

0 

-1.6200 

.7232 

i 2.6130 : 

I -2.3432 1 

-.4464 ' 


1 

1 

f 1 

- 1 . 84-77 

1 .7653 

1 2.6130 

j - 2 .ol 2 o| 

-.4142 

i 

1 2 

-3.0826 

.9176 

1 2.6130 

i - 4.0002 1 

; -.2977 

1 

1 

1 

' 3 
1 

00 

1.3065 

! 2.6130 

j ^ ! 

! 0 

4 

h 

-1.0000 

1.0000 

i 

1 CO 

-2.0000 

! -1.0000 

i 


1 1 

- 1 . 0-324 

1.0824 1 " 

- 2 .i 64 o 

j - 1.0000 



- 1.4142 

1 . 4 i 42 

! ^ 

-2 .8234 

1 - 1.0000 


I3 

-2.6130 

2 .6130 

1 CO 

-5.2260 

i - 1.0000 
» 
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TABLE CII.- INTEPPOIATION FUTTCTION 
fn,u AS CA-LOULiiTED WITH e = M = 7 


n 

i 

1 ¥ 
In.p. 

1 

0 

2.613 


1 ^ 

i -1.4l4 

i 

2 

-1.531 

I 

i 

3 

.414 

t 



i 2 

0 

-1.4i4 • 

1 

1 1 

2.696 

j 

2| 

-l.COO 

1 

3 i 
1 

-AL.531 1 

- -i 

3 

0 1 

1.082 


1 ! 

-2.4i4 

i 

1 

2 ! 

3.696 

! 

3 i 

- -i 

-1.4l4 

k ! 

0 ! 

- .500 j 

! 

1 ! 

1.082 1 

i 

» 

1 

2 j 

-l.4l4 1 


3 1 

2.613 


TA3M: cm.- VAIL7IS OF 3v n TABLE CIV.- V/iUZS OF ’o^^v 

AS CALCUL/.TSD V7IIH n = 7 AS CALCUL/iTDID ’/ITN la = 7 


i V 

1 

n ; 

; 23 v;a 

I 

I 1 

1 i 

: 0 

i 

2 i 

1 3.8284 

j 

3 i 

i 0 

i 

1 

h 1 

i .2928 

1 2 

1 

1 i 

f ■ 

i 2.0720 

1 

i 

2 i 

0 

1 

1 

3 ^ 

2.3888 



I 0 

3 

1 

0 


2 i 

1.3284 


3 i 

1 0 


4 1 

1.7022 

4 

" 1 
1 

1 1 

i .2242 

i 

2 i 

i 0 

1 

3 i 

i 3 . 15 i ^8 

j 


1 ^ 


: V : n 

21 v,V 

1 ! 1 
i 21 

i 3 ; 
; 

10.4524 

i 

2 1 1 
! 2! 

i 3! 

i ^ 

1 5.6568 

[ 

1 

i ^. ! 1 

1 ' i 2 

1 j 3 

! ^ 

4.3296 

4 1 1 
! 2 

1 ? 

1 

1 __ 

1 

! 4.0000 

i__ 
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TABLE I.- COIvIPAJRISON OF THIOKETIOAL AIJD EXPERH'DBWTAL V/vIiUES 
OF LIFT-COEVE SLOP.E AHD SPAmaSE CZMTEE OF PEESSUPE 


Wing 

parameters | 

Lift-curve slope dCL/da (deg" 

i 1 

A I 

j 

X ! 

Fallmer j 

Mutterperl i 

Weissinger | 

Experiment j 

1 ! 

1 -^ 5.2 i 

2.99 ‘ 

0.376: 

0.0419 

i 0,04005 ! 

: 0.0450 1 

0,0422 

1 -29.6 i 

4.45 I 

.405! 

.0573 1 

.0456 i 

: .0585 ! 

.0580 1 

.9 j 

4 . 47 ; 

. 542 ; 

.0633 

^ .0632 ! 

.0640 ' 

.0660 1 

31.0 j 

4.66 ! 

.442 i 

.0638 

.0615 j 

i .0631 ; 

. 0668 i 

46 A i 

3.45 •, 

i 

.418 

.0509 

.0495 i 

.0470 

' 1 

; .0538 j 

i V/ing parameters 

Spanwise center 

of pres sure 3 'Hnn ! 

I- - ■ 

! A 

AE i 

X 

Falkner 

1 

Mutterperl 

Weissinger 

Experiment | 

! -45.2 

I 

2.99 1 

0.376 

0.398 

0.385 

0.399 

0.401 j 

1 -29.6 ' 

! 4.45 1 

.405 

.408 

.362 

1 .403 

.420 i 

! .9 I 

4.47 ! 

.542 

.429 

.426 

i .425 

.433 

31.0 

4.66 I 

.442 

.439 

.434 

! .440 

1 .444 

46.4 

3.45 ! 

.418 

.446 

.438 

1 .442 

j .450 
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FIGURE I.- THE MANNER OF CONCENTRATING THE 
VORTICITY FOR THE METHODS OF FALKNER, 
MUTTERPERL AND WEISSINGER. 
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2 .4 .6 .8 /.O 

5PANWISL STATION, H 


c.) 0.9“ SWEEP 
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b.)-29.6“ SWEEP 


EXPERIMENT 

FALKNER 

MUTTERPERL 

WEISSIN6ER 



SPANWISE STATION, H 
d.)3I.O“ SWEEP 



SPANWISE STATION, f) 
e.)46;4“ SWEEP 


FIGURE 2 - DISTRIBUTION OF LOADING COEFFICIENT, ^^‘^'\/clCav ALONG THE SEMISPAN. 


LOCAL LIFT COEFnClEINT, LOCAL LIFT COEFFICIENT, LOCAL LIFT COEFFICIENT, 


NACA TN No. 1476 


77 






EXPERIMENT 

FALKNER 

MUTTERPERL 

WEISSINGER 


c.) 0.9** SWEEP 
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VIA NATIONAL ADVISOR Y 

d.)3I.O** SWEEP e.)46.4 SWEEP coMNirra rca «ik)ii*utic$ 


FIGURE 3.- DISTRIBUTION OF LOCAL LIFT COEFFICIENT, ^Vci_. ALONG THE SEMISPAN. 


